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A mi Madre 



Abstract 



In this thesis we study two important phenomenological issues in the 
context of supersymmetric gauge theories. In the Minimal Supersym- 
metric Standard Model (MSSM) we analyze properties of the neutral 
Higgs boson decays into two neutralinos, taking into account the effect 
of new quantum corrections. In the second part of our work we study in 
detail the Proton decay in the Minimal Supersymmetric SU (5) Theory. 

The main new results obtained in our studies are: 

• We compute one-loop corrections to the neutralino couplings to 
the Higgs bosons, taking into account contributions with fermions 
and sfermions inside the loop. Our analytical results are valid for 
arbitrary momenta and general sfermion mixings. For the neu- 
tralino couplings to Higgs bosons, we find in all cases corrections 
of up to a factor of two for reasonable values of the input param- 
eters. 

• The contribution of Bino-like lightest supersymmetric particles 
to the invisible decay width of the lightest MSSM Higgs boson 
might be measurable at future high-energy and high-luminosity 
e + e _ colliders, when the new quantum corrections are present. 

• The masses (Mr) of the heavy triplets T and T responsible for d = 
5 proton decay are computed, when we allow for arbitrary trilinear 
coupling of the heavy fields in E and use higher dimensional terms 
as a possible source of their masses. In this case Mr may go up 
naturally by a factor of thirty, which would increase the proton 
lifetime by a factor of 10 3 . 

• The relation between fermion and/or sfermion masses, and pro- 
ton decay is studied in detail. We find the conditions needed to 
suppress the d = 5 contributions to the decay of the proton, in 
the context of the Minimal Supersymmetric SU (5) model. 

• We point out that the Minimal Supersymmetric Grand Unified 
Theory SU (5) is not ruled out as claimed before. 
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Chapter 1 
Introduction 



In Nature the fundamental interactions are described by gauge theo- 
ries. The Standard Model based on the gauge group SU (3)c x SU{2) l x 
U(1)y explains all the properties of the electroweak and strong inter- 
actions, while Einstein's Theory of General Relativity describes the 
gravitational interaction. 

Grand Unified Theories are the main theories beyond the Standard 
Model. They explain the quantization of the electric charge, predict the 
weak mixing angle, the decay of the Proton, the bottom-tau Yukawa 
coupling unification and the existence of magnetic monopoles. At the 
same time, they provide a natural framework for understanding Baryo- 
genesis and/or Leptogenesis, and for the implementation of the see-saw 
mechanism of neutrino masses. 

It has been shown that Super symmetry (SUSY) [1], a symmetry be- 
tween fermions and bosons, plays an important role in the develop- 
ment of Unified Theories. There are many motivations for considering 
SUSY, the most important one is the possibility to cancel quadratic di- 
vergencies in the self energy of the Standard Model Higgs boson. If we 
consider radiative corrections to the Higgs mass, we see that these are 
proportional to the fundamental scale (M P i anck ~ 10 18 GeV) square, 
therefore these corrections can change its value by many orders of mag- 
nitude. This is the so-called Hierarchy Problem. Also it is possible to 
unify at the high scale M GUT ~ 10 16 GeV all the gauge coupling con- 
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stants of the Minimal Supersymmetric Standard Model [2, 3, 4, 5]. The 
possibility to break the electroweak SU(2) L x U(l)y symmetry of the 
Standard Model radiatively [6] is widely regarded as one of the main 
arguments in favor of SUSY, since it offers a dynamical explanation 
for the mysterious negative mass square of the Higgs boson. SUSY in 
the context of the Minimal Supersymmetric Standard Model provides 
us a candidate to describe the Non-Baryonic Dark Matter present in 
our Universe [7, 8, 9]. Another important motivation is the possibility 
to cancel the Tachyonic states in String Theory [10], the most popular 
scenario where all the fundamental interactions are unified. 

Another popular scenario for the solution of the Hierarchy Problem 
is large extra dimensions, which for two such new ones may be as large 
as a fraction of a mm [11, 12, 13]. In this case the field-theory cutoff 
(A F ) must be low and experiments demand: > (10 — 100) TeV. 
Clearly, one then must fine-tune (somewhat) the Higgs mass, since 



where mo and y t are the tree level Higgs mass and the top Yukawa 
coupling respectively. 

We believe this is acceptable; compared to the fine-tuning problem 
when Ap is pushed to Mpi anc k (or Mqut), this is negligible. What is 
missing in this program is some serious physical reason to have Ap so 
low. In low-energy supersymmetry, where Ap gets traded for Asusy 
(here defined as the mass difference between particles and superparti- 
cles of the MSSM). This can be as low as a few hundred GeV, therefore 
no fine-tuning whatsoever is needed. 

In our work we study two important phenomenological issues in the 
context of supersymmetric gauge theories. In the Minimal Supersym- 
metric Standard Model, we study the invisible decays of neutral Higgs 
bosons into two neutralinos at one-loop level. Proton decay in the con- 
text of Minimal Supersymmetric SU (5) is our second major objective. 




A 2 F « (few 100 GeV) 2 



(1.1) 



In the first part of the thesis we investigate the invisible Higgs decays 
into two neutralinos in the Minimal Supersymmetric Standard Model, 
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taking into account new one-loop corrections to the neutral Higgs bo- 
son couplings with neutralinos. Since the CP-odd Higgs boson does 
not couple to identical sfermions we expect that these corrections will 
be suppressed, while for the CP-even states there is the possibility to 
get large corrections. The possible impact of the corrections on the 
invisible width of the lightest CP-even Higgs boson is very important, 
because these decays could be enhanced to a level that should be easily 
measurable at future high-energy e + e~ colliders. 

In the second part of the thesis we focus on the Proton decay in the 
Minimal Supersymmetric Grand Unified Theory SU(5). We will study 
in detail the d — 5 operators contributing to the decay of the proton, 
writing down the possible contributions for each decay channel 1 . We 
point out the major sources of uncertainties in estimating the proton 
decay lifetime. We compute the masses of the color octet and weak 
triplet supermultiplets in the adjoint Higgs, in a general model where 
non-renormalizable operators are present in order to correct the relation 
between fermion masses. We study the effect of the mixings between 
fermion and sfermions in proton decay. Finally we will see if it is pos- 
sible to satisfy the experimental bounds on proton decay. 

The present thesis contains seven chapters. In the second chapter we re- 
view all the basics for Supersymmetry, we define the SUSY algebra and 
introduce all the needed tools to write down the supersymmetric ver- 
sion of gauge field theories. In chapter 3, the minimal supersymmetric 
extension of the Standard Model is introduced, all the interactions and 
relevant mass matrices for our analysis are studied. In the fourth chap- 
ter we start with the study of our first objective, the invisible Higgs 
decays into two neutralinos. We show how to compute the one-loop 
corrections, and give several numerical examples to show the effect of 
our quantum corrections. In chapter 5 we outline all the important 
aspects of the Minimal Supersymmetric Grand Unified Theory SU(5). 
In Chapter 6, we study our second important phenomenological issue, 
Proton decay. We discuss all the relevant operators contributing to 



1 Here d refers to the mass dimension of the operator, not to the dimension of 
spacetime. 
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the decay of the proton in supersymmetric theories, and we focus our 
analysis in the Minimal Supersymmetric SU(5). Finally in Chapter 7 
we conclude, pointing out possible future directions. 



Chapter 2 

Basics of Supersymmetry 



2.1 SUSY Algebra 

Supersymmetry is a symmetry between fermions and bosons, which is 
generated by a fermionic generator Q. 

Q 

^ fermionic *>' ^bosonic 

In general we could define a Supersymmetric Field Theory, as a theory 
which is invariant under SUSY transformation: 

S S usy(Q)S = Ssusy(Q) J d D xC(^!) = 

With the usual Poincare and internal symmetry algebra, it is possible 
to define the Super-Poincare Lie algebra, which contains the additional 
SUSY generators Q l a and Qj., where Q\ = (Q'JiU] [15]: 

[^,^]=0 (2.1) 
[Pp,M pa \=i(g w P a -g lia P p ) (2.2) 

[M^, M pc \ = i(g up M^ - g va M w - g w M va + g^M vp ) (2.3) 

5 



6 



CHAPTER 2. BASICS OF SUPERSYMMETRY 



[B r ,B a ]=i<* a B t (2.4 
[B r ,P,} = (2.5 
[B r ,M fia ] = (2.6 
[Qa,P,] = [Ql,P,] = (2.7 
[QlM^^^aM/s (2-8 

[ft,^] = -^^V)f (2-9 

[Q i a ,B r ] = (b r ))Q> a (2.10 

[Ql,B r ] = -Qi(b r )) (2.11 

{Qj,,gj} = 2^K)^ (2.12 

{QLQj} = 2e^' (2.13 

Z l3 = a^Br, Z'i = 4 (2.14 

{0:,.0\\ 2< ni Z' ! (2.15 

[Z^, anything] = (2.16 



where a , a — 1, 2; i,j — l,2,...,N, with iV as the number of super- 
symmetries. 

Here P M is the four-momentum operator, My and M j are the angular 
momentum and boost operators respectively, B r the internal symmetry 
generators, g a p is the metric, c l rs and a\^ are structure constants and Z^ 
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are the so-called central charges; a,a,f3,f3 are spinorial indices. In the 
simplest case one has one spinor generator Q a (and the conjugated one 
Qa) that corresponds to an ordinary or N=l supersymmetry. It is has 
been proved that the Super-Poincare Lie algebra contains all possible 
symmetry generators for symmetries of the S-matrix. It is the so-called 
the Coleman-Mandula Theorem [16]. 

There are many important conclusions coming from the SUSY alge- 
bra. We see from equations 2.8 and 2.9, that the SUSY generators 
change the spin by a half-odd amount and change the statistics. While 
from equation 2.7 we can conclude that a fermionic (or bosonic) field 
and its superpartner in a theory with exact supersymmetry must have 
the same mass. It is the reason why SUSY must be broken in order to 
get a realistic spectrum in particle physics. 

2.2 Superspace and Superfields 

An elegant formulation of supersymmetric transformations and invari- 
ants can be achieved in the framework of superspace [17]. Superspace 
differs from the ordinary Euclidean (Minkowski) space by adding two 
new coordinates, 6 a and 9 a , which are Grassmannian, i.e. anticom- 
muting, variables 

{e a ,e p } = o, {e a ,e $ } = o, e 2 a = o, e 2 .=o, a,p,a,p = 1,2. 

Thus, we go from space to superspace 

Space =>- Superspace 
x^, Qq, , 6 a 

A SUSY group element can be constructed in superspace in the same 
way as an ordinary translation in the usual space 

G(x,e,0) = e i (- x * P » + 6 Q + 8Q) (2.17) 
It leads to a supertranslation in superspace 

3>n -> Xf, + iOa^E - ieaj) 
9_ -> 9_ + e (2.18) 
6 -> 9 + e 
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where e and e are Grassmannian transformation parameters. From 
eq.(2.18) one can easily obtain the representation for the supercharges 
acting on the superspace 

Qa = -^-i<rL0% Q* = -g^ + i0aPLd ll (2.19) 

We are now ready to introduce the superfields. The superfields can 
be defined as functions in Superspace, F = F(x,6,6). However, these 
superfields are in general reducible representations of the SUSY algebra. 
To get an irreducible one, we define a chiral superfield $ which obeys 
the equation: 

d 

£>$ = where D = — -= — ida^d, (2.20) 

is a superspace covariant derivative. For the chiral superfield, the Grass- 
mannian Taylor expansion looks like (y = x + idaO) 

$(y,e) = A(y) + V26^(y) + 00F(y) 

= a(x) + i9a^A(x) + -eeemA(x) 

+ V2ej/j(x)--^=eed^(x)a^e + eeF(x) (2.21) 

V2 

The coefficients are ordinary functions of x, being the usual fields. 
There are two physical fields, a bosonic one A and a fermionic ip, while 
F(x) is an auxiliary field without physical meaning, needed to close the 
SUSY algebra. 

Under SUSY transformation the fields convert into one another 

5 £ A = V2ei> 

5 £ ip = iV2a^ed^A + V2eF (2.22) 
5 £ F = iV^ea^d^ 

Note that the variation of the F-component is a total derivative, i.e. , 
with appropriate boundary conditions it vanishes when integrated over 
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the space-time. 

One can also construct an antichiral superfield & obeying the equation: 

Q 

D& = 0, with D = — + ia^dda 

06 

The product of chiral (antichiral) superfields $ 2 , $ 3 , etc is also a chiral 
(antichiral) superfield, while the product of chiral and antichiral ones 
$t$ is a general superfield, it is not a chiral superfield, its 6666 com- 
ponent transforms under SUSY as a total divergence. 
To construct the gauge invariant interactions, one needs a real vector 
superfield, which is defined as V = V\ The explicit form of V is: 



V(x,6,6) = C(x)+i6 X (x) -i6x(x) 

+ -66[M(x) + iN(x)\ - -66[M(x) - iN(x)} 

Z z 

- da^dv^x) + i666[\{x) + l -^d^x(x)] 

- iBM[\ + \<J^x{x)\ + \6666[D{x) + \uC{x)\ 

(2.23) 

The physical degrees of freedom corresponding to a real vector super- 
field V are the vector gauge field and the Majorana spinor field A. 
All other components are unphysical and can be eliminated. 

Under the Abelian (super)gauge transformation the superfield V is 
transformed as V — > V + $ + where $ and $' are some 
chiral superfields. In components it looks like [14] 

C -> C + A + A* 

X -> X - iv^V, 

M + iN -> M + iN -2iF 

v, - — idn(A — A*) (2.24) 

A — > A 

D — > £> 
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where A* is the complex conjugate of A. According to eq.(2.24), one 
can choose a gauge (the Wess-Zumino gauge) where C = x = M = 
N = 0, leaving one with only physical degrees of freedom except for 
the auxiliary field D. In this gauge 

v = -e^ev^x) + i999\{x) - mo\(x) + -9999D(x) 

V 3 = etc. (2.25) 

2.3 Supersymmetric Lagrangians 

Using the rules of Grassmannian integration: 

J de a = o Je a de, = 5 a , 

we can define the general form of a SUSY and gauge invariant la- 
grangian as [14]: 

£ Y susy = \jd 2 9Tr(W a W a ) + ^Jd 2 9Tr(W a W a ) (2.26) 
+ J d 2 9d 2 9 $J a {e 9V ) a b <&\ + J d 2 9 W{$i) + J d 2 9 W($,) 

$i are chiral superfields which transform as: 

$. _^ e - i9A <s>i 

and 

where, both A and V are matrices: 

Aij = r«A a Vn = r*V a 

with r a the gauge generators. The supersymmetric field strength W a 
is equal to 
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and transforms as: W — > e lA W^e 



W is the superpotential, which should be invariant under the group 
of symmetries of a particular model. 

In terms of component fields the above Lagrangian takes the form [18] 



&msY = ~\f; v F^ -aV^A« + l -D a D a 

+ (d^Ai - igv;T a A t )\d,A, - igv^r'AJ ~ ^{d^i - igv a »T a ^) 
- D a A\r a A l - iV2Alr a \ a ^ + i^^AX + F-Fi 

Integrating out the auxiliary fields D a and F i: one reproduces the usual 
Lagrangian. 

Contrary to the SM, where the scalar Higgs potential is arbitrary and is 
defined only by the requirement of the gauge invariance, in supersym- 
metric theories it is completely defined by the superpotential. It con- 
sists of the contributions from the D-terms and F-terms. The kinetic 
energy of the gauge fields yields the ^D a D a term, and the matter-gauge 
interaction yields the gD a r^A*Aj one. Together they give 

C D = l -D a D a + gD a T^A*Aj (2.28) 

The equation of motion reads 

D a = -gr^Aj (2.29) 

Substituting it back into eq.(2.28) yields the D-term part of the poten- 
tial 

C D = --D a D a Vb = -D a D a (2.30) 
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where D is given by eq.(2.29). 

The F-term contribution can be derived from the matter field self- 
interaction. For a general type superpotential W one has 

C F = F*Fi + (— + h.c.) (2.31) 

Using the equations of motion for the auxiliary field Fj 

yields 

C F = —F*Fi ^V F = F*Fi (2.33) 

where F is given by eq.(2.32). The full potential is the sum of the two 
contributions 

V = V D + V F (2.34) 

Thus, the form of the Lagrangian is constrained by symmetry require- 
ments. The only freedom is the field content, the value of the gauge 
coupling g, Yukawa couplings y^ and the masses. Because of the renor- 
malizability constraint V < A A the superpotential should be limited by 
W < $ 3 . All members of a supermultiplet have the same masses, i.e. 
bosons and fermions are degenerate in masses. This property of SUSY 
theories contradicts the phenomenology and requires supersymmetry 
breaking. 



2.4 SUSY Breaking 

Since the supersymmetric algebra leads to mass degeneracy in a super- 
multiplet, it should be broken to explain the absence of superpartners 
at accessible energies. There are several ways of supersymmetry break- 
ing. It can be broken either explicitly or spontaneously. Performing 
SUSY breaking one has to be careful not to spoil the cancellation of 
quadratic divergencies which allows one to solve the Hierarchy problem. 
This is achieved by spontaneous breaking of SUSY. 
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It is possible show that in SUSY models the energy is always nonneg- 
ative definite. According to quantum mechanics the energy is equal 
to: 

E=(0\H\0) (2.35) 
where H is the Hamiltonian and due to the SUSY algebra: 

{Q«,0 / j} = 2K)^P /l (2.36) 
taking into account that Tr(a^P fl ) = 2P one gets 

E = 7 E < 0|{Q«> 0J|0 >= j £ ||Q Q |0 > || 2 > (2.37) 

4 o=l,2 4 a 

Hence 

£ =< 0| H |0 >^ i/ and only if Q a \0 >^ 

Therefore, supersymmetry is spontaneously broken, i.e. the vacuum is 
not invariant under Q (Q a \0) ^ 0), if and only if the minimum of the 
potential is positive (i.e. E > 0) . 

Spontaneous breaking of supersymmetry is achieved in the same way 
as electroweak symmetry breaking. One introduces a field whose vac- 
uum expectation value is nonzero and breaks the symmetry. However, 
due to the special character of SUSY, this should be a superfield whose 
auxiliary F or D component acquires nonzero v.e.v.'s. Thus, among 
possible spontaneous SUSY breaking mechanisms one distinguishes the 
F and D ones. 

i) Fayet-Iliopoulos (D-term) mechanism [18]. 

In this case the, the linear D-term is added to the Lagrangian 



A£ = ty\oe-ee = i j d 2 6 d 2 9 V 



(2.38) 
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It is U(l) gauge and SUSY invariant by itself; however, it may lead to 
spontaneous breaking of both of them depending on the value of £. The 
drawback of this mechanism is the necessity of U{1) gauge invariance. 
It can be used in SUSY generalizations of the SM but not in GUTs. 
The mass spectrum also causes some troubles since the following sum 
rule is always valid 

STrM 2 = ]T(-1) 2J (2J + l)m 2 j = (2.39) 
j 

which is bad for phenomenology. 

ii) O'Raifeartaigh (F-term) mechanism [18]. 

In this case, several chiral fields are needed and the superpotential 
should be chosen in such way that trivial zero v.e.v.s for the auxiliary 
F-fields are forbidden. For instance, choosing the superpotential to be: 

W($) = A$ 3 + m$i$ 2 + g$ 3 ®l (2.40) 
one gets the equations for the auxiliary fields 

F* = mA 2 + 2gA 1 A 3 (2.41) 
F; = mA 1 (2.42) 
F* = X + gAl (2.43) 

which have no solutions with (Fj) = and SUSY is spontaneously bro- 
ken. 

The drawback of this mechanism is, that there is a lot of arbitrari- 
ness in the choice of potential. The sum rule (2.39) is also valid here. 

Unfortunately, none of these mechanisms explicitly works in SUSY 
generalizations of the SM. None of the fields of the SM can develop 
nonzero v.e.v.s for their F or D components without breaking SU(3) 
or U(l) gauge invariance since they are not singlets with respect to 
these groups. This requires the presence of extra sources for sponta- 
neous SUSY breaking [19, 20, 21, 22, 23, 24]. 



Chapter 3 
The MSSM 



The Standard Model (SM) describes with a very good precision all 
electroweak and strong processes. It is based on gauge invariance under 
the symmetry group: 

G SM = SU(3) C x SU(2) L x U(1) Y (3.1) 

and its partial spontaneous symmetry breaking. In Table 1 we show 
all its constituents, the elementary fermions (quarks and leptons), the 
scalar Higgs boson, the gauge bosons, and their transformation prop- 
erties under Gsm- We use the relation Q = T 3 + Y/2, where Q, T 3 and 
Y are the electric charge, isospin and hypercharge respectively. 

The lagrangian of the Standard Model has the following form: 

£-SM £ fermions t~-gauge bosons ~^ ^scalars ^Yukawa 

The explicit form of the SM lagrangian is well known, for our objectives 
we will write explicitly only the expression of the Yukawa interactions: 



C Y uka W a = (U l d% h% $ 4+(^ d l ) L hfy <7 2 $* 4+(P j e l ) L h% $ 4+/I.C 

(3.3) 

u\ d\ h u and h d are the quarks with isospin 1/2 and —1/2, and their 
Yukawa matrices respectively, e % and h e stand for the charged leptons 
and their Yukawa matrices, while v l are the neutrinos for each family. 
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The subscripts L and R refer to right and left chirality respectively, 
while i and j are the generation indices. 



Table 1. Standard Model Particles 



Quarks: 



u 
d a 



R 



where a — 1, 2, 3 (colors) 
Leptons: 



eR 
Scalars: 



Gauge bosons: 
Gl with a = 1,2. .8 
withb= 1,2,3 



6° 



S77(3) c , 5C/(2) L ,C/(l) y 

(3c,2 L ,l/3) 

(3c, U, 4/3) 
(3 C , 1 L , -2/3) 



(1c,2 l ,-1) 
(lc,U,-2) 



;ic,2 L ,l) 



(8c,U,0) 
(lc,3 L ,0) 
(1c,U,0) 
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Note that we can write these terms using only one scalar field $, which 
after electroweak symmetry breaking can generate mass for all the 
quarks and leptons in the Standard Model. 

An important free parameter in the Standard Model is the Weinberg 
angle 9 W , which is defined as: 



where gu(i) an d gsu{i) are the gauge couplings for the U(l) and SU{2) 
gauge groups. 

As was mentioned above, the standard model has an extremely econom- 
ical Higgs sector, which accounts for all the particle masses. Baryon 
(B) and Lepton (L) numbers are automatically conserved and it is an 
anomaly free Quantum Field Theory. However not all is perfect, at 
present there is no evidence for Higgs. 

A dramatic problem in the Standard Model is present in the Higgs 
sector. If we consider radiative corrections to the Higgs mass, we see 
that it has quadratic divergencies, which can change its value by many 
orders of magnitude [25]. In Fig.3.1 we show two of the problematic con- 
tributions, these are the contributions with fermions and gauge bosons 
inside the loops. The only known way to cancel these divergencies is 
supersymmetry. SUSY automatically cancels quadratic corrections to 
all orders of perturbation theory. This is due to the contributions of 
superpartners of ordinary particles. The contribution from boson loops 
cancel those from the fermion ones because of an additional factor (-1) 
coming from Fermi statistics. 

The first line of Figure 3.1 shows the contribution of an SM fermion and 
its superpartner. The strength of interaction is given by the Yukawa 
coupling A. The second line represents the gauge interaction propor- 
tional to the gauge coupling constant g with the contribution from the 
gauge boson and its superpartner. 



sin 6 



9u(i) 



(3.4) 



'w = 
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Figure 3.1: Cancellation of quadratic divergences (from reference [26]) 

In both cases, cancellation of quadratic terms takes place. This cancel- 
lation occurs up to the SUSY breaking scale, A SUSY , since 

E ™ 2 " E m " = a2 susy (3-5) 

bosons fermions 

which should be around ~ 1 TeV to make the fine-tuning natural. 
Indeed, let us take the Higgs boson mass. Requiring for consistency of 
perturbation theory that the radiative corrections to the Higgs boson 
mass do not exceed the mass itself gives 



l07T z 



2., d « A2 _ . n,r2 



So, if M h ~ 10 2 GeV one needs A S usy ~ 10 3 GeV in order that the 
relation (3.6) is valid. Thus, we again get more or less the same rough 
estimate of A SUSY ~ 1 TeV as from the gauge coupling unification. 
Two requirements match together. However as we mentioned before, 
Asusy could be in the range 1 — 10 TeV if we accept some fine-tuning. 
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19 



3.1 Particles and their Superpartners 

In the MSSM, we add a superpartner for each SM particle in the same 
representation of the gauge group. Usually we use the SUSY operators 
in the left-chiral representation, therefore it is convenient to rewrite 
the SM particles (Table 1.) in the left-chiral representation and define 
the superpartners accordingly. This leads to the superfield formalism, 
which makes it easier to construct SUSY invariant lagrangians. In this 
case we have to introduce for each SM particle one superfield, which 
contains the SM particle, its superpartner and an auxiliary unphysical 
field. In Table 3 we show the third generation of the SM particles, their 
superpartners and the superfields needed to write our lagrangian. Note 
that we have an extended Higgs sector and the color index is omitted. 
In Table 2 we show the names of the superpartners. 



Table 2. Names of superpartners. 



Matter Fermions 



Sfermions 



(quarks, leptons) 
s = 1/2 
Gauge Bosons 
(W ± , Z, 7, gluons) 
s = 1 
Higgs 
s = 



(squarks, sleptons) 
s = 
Gauginos 
(Wino, Zino, photino, gluinos) 
s = 1/2 
Higgsinos 
s = 1/2 



For example the superpartner of the top quark is called stop, of the 
photon the superpartner is the photino, and similarly for the other 
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particles. However, in the Higgs sector we must add another new Higgs 
boson and its superpartner to write the Yukawa interactions needed to 
generate masses for all quarks and to obtain an anomaly free model. 



3.2 MSSM Lagrangian 

We can divide the lagrangian of the minimal extension of the Standard 
Model into two fundamental parts, the SUSY invariant and the Soft 
breaking term [27]: 

£>MSSM = C-susy + £soft (3.7) 
In general we can write the SUSY invariant term as: 

£-SUSY = £gauge + Qeptons + ^quarks + ^Higgs + J d 2 W + h.C (3.8) 

Defining the content of the MSSM in Table 3, we can write the different 
terms of the lagrangian. The term C gauge has the following form: 

Cgauge = \j d 2 0[2 Tr(W 3 W 3 ) + 2 Tr(W 2 W 2 ) + W l W^ (3.9) 
with: 

Wl = ~D Dexp(-G 3 ) D a exp(G 3 ) (3.10) 

Gs = j:^G a 3 (3.11) 

a=l Z 

W 2 = ~D Dexp(-G 2 ) D a exp(G 2 ) (3.12) 
<?2 = EirG$ (3-13) 

6=1 Z 
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Table 3. Content of the MSSM. 
Superfields 



Vector Superfields 



Bosonic Fields Fermionic Fields Gsm 



G3 with a — 1, 2. .8 G a (8 C ,U,0) 

G\ W* with b= 1,2,3 (lc,3 L ,0) 
Gi ^ B (1c,1l,0) 

Chiral Superfields 
Leptons 



.By \ey \e. 

£ C e c e c (lo,U,2) 

Quarks 

g=(£) 3= (|) (3) (3c,2„l/3) 

f/ C S c « c (3c, 1 L , -4/3) 

B c d c rf c (3 c ,1l,2/3) 

Higgs 

MCO (I) {lc ^- 1] 

H H 2 =( H J ) (f ) (1 C ,2 L ,1) 
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and 

W l a = --^DDD a Gi (3.14) 

where A a and o a are the Gell-Mann and Pauli matrices respectively, 
and G's are defined in Table 3. 

The SUSY covariant derivatives D and D are used in the left chiral 
representation: 

Dl = ^ + 2ia»9d, (3.15) 

and 

D L = -| (3.16) 

The lagrangians for gauge interactions of leptons, quarks and the Higgs 
bosons are: 

Cleans = J d 2 9 d 2 9 Lt exp(2g 2 G 2 + gJ^Gi) L + 

+ J d 2 9 d 2 9 E c] ex^{2g 2 G 2 + gJ^G^) E c (3.17) 
C quarks = J d 2 6d 2 9 Q f exp(2g 3 G 3 + 2g 2 G 2 + gJ^-G 1 ) Q + 
+ J d 2 9 d 2 9 U ct exp^^pd - g 3 (\ a TG a 3 ) U c + 
+ J d 2 9 d 2 9 D c] exp(^ 1 ^G 1 - g 3 {\ a YG a 3 ) D° (3.18) 
C mggs = J d 2 9 d 2 9 & exp(2g 2 G 2 + gj^d) H + 

+ f d 2 9 d 2 9 H ] exp(2g 2 G 2 + ffi^d) # (3.19) 

where g 3 , g 2 , gi are the SU(3), SU(2) and U(l) coupling constants, re- 
spectively. The Yi represent the hypercharges of the different super- 
fields. 
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We can write the superpotential as the sum of two terms, W=Wr+Wnr- 
The first conserves lepton (L) and baryon (B) numbers: 



W R = t tJ [-^H l H ] + H l VY E E c + H l Q ] Y D D c + iPQ%£/ c ] (3.20) 



where eij is the antisymmetric tensor, \i the Higgs mass parameter and 
Yu, Y D and Y E are the different Yukawa matrices. The term Wnr, 
which explicitly breaks L and B numbers, is: 



W NR = t tj [-^H l L j + \VVE C + \'VQ j D c ] + \"D C D C U C (3.21) 



In Wnr the first three terms break lepton number, while the last term 
breaks baryon number. From these terms we find d — 4 operators 
contributing to the decay of the proton, which will be analyzed in the 
next chapters. In the last section of this chapter we will analyze the 
R-symmetry related with the L and B number conservation and its im- 
plications. 

SUSY is broken explicitly if we introduce the following terms: 



-C ao ft = m 2 1 \H 1 \ 2 + m\\H 2 \ 2 + m\ 2 {E x E 2 + H{H* 2 ) 

+ Q ] M 2 ~Q + u^m\ c u c + dP^m 2 ~ c d c + ZtMfZ + m\ c E 
+ [H 2 Q{YuAu)u c + H^YdAd)! + H X L{Y E A E )E C + h.c] 



Note that in order to describe SUSY breaking we introduce many free 
parameters, and several terms have mass dimension less than 4 (super- 
renormalizable, but not SUSY invariant). The different mass terms 
remove the degeneracy between particles and their superpartners. 




(3.22) 



3.3 Neutralinos and Charginos 



Once SU(2) L x U(1) Y is broken in the MSSM, fields with different 
SU(2)l x U(1)y quantum numbers can mix, if they have the same 
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SU(3)c x U(l) em quantum numbers, and the same spin. 

The neutralinos are mixtures of the B, the neutral W 3 and the two 
neutral Higgsinos. 

The mass term for the neutralinos is equal to: 

-CJ SS = U^f M N ^ + h.c (3.23) 

If we define the physical states as Xi — Nijip®, the diagonal mass matrix 
is M D = N*M N Nl 

In general these states form four distinct Majorana fermions, which 
are eigenstates of the symmetric mass matrix [in the basis (B, W s , 
flf)][27]: 



M 



N 



(3.24) 



Mi -m z s w cp m z s w sp 

M 2 m z c w cp -m z c w S(3 

—m z swC/3 m z cwCp — // 

m z swsp —m z cwSj3 —fi 

where Mi and M 2 are the SUSY breaking masses for the U(l)y and 
SU(2)x, gauginos, ji is the higgsino mass parameter, and sp = sin/5, 
s w = sm9 w , etc. 

We will assume that all soft breaking parameters as well as fx are real, 
i.e. conserve CP. We can then work with a real, orthogonal neutralino 
mixing matrix iV if we allow the eigenvalues m^o to be negative. 

This matrix can be diagonalized analytically, but the expressions of 
the neutralino masses and the iV^ matrix elements are rather involved. 
However, if the entries in the off-diagonal 2x2 submatrices are small 
compared to the diagonal entries, one can expand the eigenvalues in 
powers of m z [28] : 



m 



M i " -r-^TH ( M i + /^) s w (3-25) 
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A2 



A3 



X4 



~ M 2 - 



(M 2 + lls w ) <? w 



rrtr 



-fi 



Ll + 



m 2 z (l + sip) 



+ 



-w 



Ll + M 1 Ll + M 2 



s 2 



+ 



-w 



Li- Mi Li - M 2 



(3.26) 
(3.27) 
(3.28) 



In our analysis, we are interested in the situation > Mi, M 2 and 
Li 2 ^> m|. In this case the lighter of the two neutralinos will be gaugino- 
like. If I Mi I < |M 2 |, the lightest state will be bino-like, and the next- 
to-lightest state will be wino-like. The two heaviest states will be 
dominated by their higgsino components. The components of the mass 
eigenvectors can also be expanded in powers of m z . We find for the 
bino-like state: 



N n 
N13 
N n 

N n 



1 + (N 12 /N u ) 2 + (iV 13 /iVn) 2 + (N lA /Nn? 



-1/2 



m? z s w c w (s 2/3 /i + Mi 
fi' 

m z s w 



-m z s w 



Ml (Mi - M 

{SfjLL + CpMx) 
Li 2 - Ml 
(cpfi + spMx) 



+ 0{m\) 



fi 2 - Ml 



+ 0{m 2 z ) 
+ 0{m\) 



The corresponding expressions for the wino-like state read: 



N 22 
^21 

^23 
A^22 
^4 
N 22 



1 + (iV 2 i/iV 22 ) 2 + (N 2 ,/N 22 ) 2 + (N 24 /N 22 f 



-1/2 



m 2 z s w c w (s 2f} Li + M 2 ) 



Mi (M 



-m z c w - 



2 -Mi) 

(spLi + CpM 2 ) 



M 2 



= m z c w 



fi- - ivi 2 

fi 2 - M| 



+ 0(m%) 
+ 0{m\) 



+ 0(m%) 



(3.29) 
(3.30) 

(3.31) 

(3.32) 

(3.33) 
(3.34) 

(3.35) 

(3.36) 



Note that the higgsino components of the gaugino-like states start at 
0(mz), whereas the masses of these states deviate from their — > 00 
limit (Mi and M 2 ) only at 0(m 2 z ). 
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The charginos are mixtures of the W ± and H ± . The chargino mass 
matrix [in the basis (W ± ,H i ) ] [27]is: 



M 



c 



M 2 \f2m w sp 
y/2m w C(3 pL 



(3.37) 



if we expand in powers of mw, the two chargino masses are: 



m~± ~ M 2 

Al 



m 



w 



(M 2 + ns 2p ) 



(3.38) 



m ^ ~ * + ^ - M| (M2 ' 2/3 + ^ (3 ' 39) 

so that for \fi\ — > oo, the lightest chargino corresponds to a pure wino 
state while the heavier chargino corresponds to a pure higgsino state. 

Usually the neutralino is considered the lightest supersymmetric parti- 
cle (LSP) in models where R-parity is conserved. It has been realized 
many years ago that they are good candidates to describe the Non- 
Baryonic Dark Matter present in the Universe [7, 8, 9]. There is no 
direct experimental limit for the neutralino mass, however from LEP 
experiments we know that the chargino mass m ± must be bigger than 

103.5 GeV[29]. 



3.4 Squarks and Sleptons 

After the electroweak symmetry breaking, several terms in the MSSM 
lagrangian contribute to the sfermion mass matrices. Ignoring flavor 
mixing between sfermions, the mass matrix for charged matter sfermion 
is[in the basis (/, f c )} [27]: 



M 2 = ( mj + m 2 LL m f A f \ 
f V m f m) + m\ R J 



(3.40) 
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with 

m\ L = m 2 j + — Qf sin 2 9 W ) m| cos 2(3 
m 2 RR = m 2 j c + Q f sin 2 9 w m 2 z cos 2(3 (3.41) 
A f = A f -/i(tan/5)- 2/ 3 

where / represents the different charged fermions u\ d l and e\ 

The charged sfermions mass matrices are diagonalized by 2 x 2 rotation 
matrices described by the angles 9j, which turn the current eigenstates, 

/ and f c , into the mass eigenstates fi and f2] the mixing angle and 
sfermion masses are then given by 

,2 



2 2,1 

m- f = nit H — 

/l,2 / 2 



sm2e f = 2 2 mfAf 2 , cos2g f -= m ^ m f , (3.42) 

(3.43) 



m ; — m'- J ra^ — mi, 

71 72 /l 72 



2l + m RR T \J(m 2 LL - m 2 RR ) 2 + Am 2 A 2 



The physical states are defined as: 

/i = / cos flj + / c sin 9 j (3.44) 



and 

/ 2 = -/sin flj + / c cos 9 j (3.45) 
In the case of sneutrinos we have: 

M 2 = Ml + -m| cos 2(3 (3.46) 

Note the contributions of the different soft breaking parameters in the 
mass matrices. 
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3.5 Higgs Bosons 

As we have mentioned before, the existence of a scalar Higgs boson is 
the main motivation to introduce SUSY. In the MSSM the tree-level 
Higgs potential is given by [27] [30] : 



V Higg s = ™ffj#i| 2 + m 2 H2 \H 2 \ 2 + m\ 2 {H x H 2 + h.c)+ (3.47) 



2 2 2 

+ 9i±92q Hi \2 _ |# 2 |2 )2 + ?l\ H l H2 f (3.4 8 ) 
8 2 

where m 2 H . = mf + |/i| 2 with i — 1,2 

Note that in this equation the strength of the quartic interactions is 
determined by the gauge couplings. 

After electroweak symmetry breaking, three of the eight degrees of free- 
dom contained in the two Higgs boson doublets get eaten by the W ± 
and Z gauge bosons. The five physical degrees of freedom that remain 
form a neutral pseudoscalar boson A , two neutral scalar Higgs bosons 
h° and H°, and two charged Higgs bosons H + and H~. 

The physical pseudoscalar Higgs boson A is a linear combination of 
the imaginary parts of H® and H 2 , which have the mass matrix [ in the 

basis (^,^#2)]: 

m 2 A0 = trMf = -2m\ 2 / sin(2/3) (3.50) 

The other neutral Higgs bosons are mixtures of the real parts of 
and H®, with tree-level mass matrix \—^ L , ~§ 2 ~] : 



-m\ 2 tan f3 + m 2 z 
m i2 ~~ k m 'z sn: 



; z cos' 
-"-^1 sin 2(3 



' (3 m\ 2 — \m 2 z $ni2(3 
—m 2 2 co * P + m l s™ 2 



(3.51) 



3.6. THE R-SYMMETRY AND ITS IMPLICATIONS 
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In this case the eigenvalues are: 
1 

2' 



m H° h° = o [ m A0 + m l =•= V ( m A» + m z) 2 ~ 4m 2 AO m 2 z cos 2 (2/3)] (3.52) 



Explicitly the mass eigenstates are: 

j G° = cos/5 im if? + sin/3 7m # 2 °], Goldstone boson -> Z°, 

| A = ^[sin/3 Jm + cos/3 Im Jf£], Pseudoscalar Higgs, 

j G + = [- cos (3 (iff ) * + sin (3 # 2 +] , Goldstone boson ->W + , 

\ H + = ^[sin/3 (iff)* + cos/3 ff 2 + ], Charged Higgs, 

{h° = -^[— sin a Re if° + cos a f?e if"], light Higgs, 

H° = ^[cosa Re ff° + sin a Re if 2 ], heavy Higgs, 



where the mixing angle a is given by: 

nn ( m 2 ,o + Ml 
tan2a = - tan 2/3 — £ -f 

Wa<> - 

From these equations we can see that at tree level, the MSSM predict 
that m h o < mz, however when we consider one- loop corrections, the 
mass of the light Higgs boson is modified significantly. For example 
assuming that the stop masses do not exceed 1 TeV, m h o < 130 GeV 
[31]. 

3.6 The R-symmetry and its Implications 

In the Standard Model the conservation of Baryon (B) and Lepton (L) 
number is automatic, this is an accidental consequence of the gauge 
group and matter content. In the MSSM, as we showed in the second 
section of this chapter, we can separate the most general gauge invari- 
ant superpotential into two fundamental parts, where the first term 
conserves B and L, while the second breaks these symmetries. 



(3.53) 



30 



CHAPTER 3. THE MSSM 



In the MSSM, B and L conservation can be related to a new discrete 
symmetry, which can be used to classify the two kinds of contributions 
to the superpotential. This symmetry is the matter-parity, defined as: 

M = (-lf B -V (3.54) 

Quark and lepton supermultiplets have M = — 1, while the Higgs and 
gauge supermultiplets have M — +1. The symmetry principle in this 
case will be that a term in the lagrangian is allowed only if the product 
of the M parities is equal to I. 

The conservation of matter-parity as defined in equation (3.54), to- 
gether with spin conservation, also implies the conservation of another 
discrete symmetry called R-parity, defined such that it will be +1 for 
the SM particles and -1 for all the sparticles: 

R=(-1) 2S M (3.55) 

These two symmetries are equivalent, since in the superpotential only 
the scalar fields get v.e.v. However only M commutes with the SUSY 
generators. 

Now if we impose the conservation of R, we will have some important 
phenomenological consequences in SUSY models: 

• The lightest particle with R — — 1, called the lightest supersym- 
metric particle (LSP), must be stable. 

• Each sparticle other than the LSP must decay into a state with 
an odd number of LSPs. 

• Sparticles can only be produced in even numbers from SM parti- 
cles. 

• There are not d = 4 operators contributing to the decay of the 
proton. 

It is important to mention that the conservation of R parity is predicted 
in a large class of Supersymmetric Grand Unified Theories as Minimal 
SUSY 50(10) [32, 33]. 



Chapter 4 

SUSY decays of neutral 
Higgs bosons 



In this chapter we will analyze different aspects related to supersym- 
metric decays of the Higgs bosons in the Minimal Supersymmetric Stan- 
dard Model (MSSM), studying the Higgs decays into two neutralinos. 
In particular we will compute and show the effect of new loop correc- 
tions to the Higgs-neutralino-neutralino couplings and to the invisible 
branching ratios. 



4.1 Higgs decays in the MSSM 

In order to provide a complete analysis of the most important aspects of 
Higgs decays in the Minimal Supersymmetric Standard Model (MSSM), 
we will start with the properties of the SM Higgs boson. The Higgs 
mass Tn\ SM = \\v 2 is a free parameter, since A is unknown at present. 
However the theory predicts the Higgs couplings to fermions and gauge 
bosons as: 

SM _ 171 7 SM _ ^ m V (a I \ 

9hff ~ — 9hvv ~ y^- L ) 

where / is used for any fermion, and V for W ± and Z°. In Fig 4.1 
the different SM Higgs branching ratios versus rrih SM is plotted. From 
this figure we can appreciate that there are two important intervals, 
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Figure 4.1: The decay branching ratios (left) and the total decay width 
(right) of the SM Higgs boson as a function of its mass, [from reference [34]] 



for rrih SM < 135 GeV the most important channel is hsu — ► b b with 
branching ratio close to ~ 90%, while for m hsM > 135 GeV the dom- 
inant decay mode is h SM — > W W* (where one of the gauge bosons 
may be virtual). This behaviour is easy to understand if we take a 
look at the couplings listed above. There are other important chan- 
nels such as hsM ^ r + r~ , hsM c c, and at one- loop we have the 
decays hsu 9 9 and h S M -> 7 7, which are important for Higgs 
searches [31]. In the low mass range, the Higgs width is very narrow, 
with T < 10 MeV, but increasing m hsM we reach 1 GeV at the Z° Z° 
threshold (see Fig 4.1). 

The same analysis in the context of the MSSM is more difficult, due to 
the presence of three neutral Higgs bosons and one charged pair. There 
are new decay modes which modify the branching ratios of all the chan- 
nels, in particular the decay into supersymmetric particles could play 
an important role [35, 36, 37, 38]. 

In order to understand how the branching ratio are modified, we list the 
couplings of the neutral Higgs bosons to / / relative to the Standard 
Model values: 
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for the light CP-even Higgs h°: 

cos/3 sin/3 



t n t t / \ sm a 7 n cos a / ^ r.\ 

/r dj dj( or J : /i Uj tij : — — — (4.2) 



for the heavy CP-even H° 

„ n , - . . cos a t sin a: 

#° d, or e t €i : H° m u % : — — 4.3 

cos p sm p 

and for the CP-odd A we have: 

A di d~i( or e,i e"j ) : 75 tan /3 A Mj Uj : 75 cot /3 (4.4) 



while the couplings of the two CP-even Higgs bosons to W ± and Z° 
pairs are given by: 

9b?vv = Qh^vv sin (/5 - a ) 9h«vv = Qh^vv cos (P - a ) ( 4 - 5 ) 



From the couplings listed above, we note that there are two new pa- 
rameters which will play an important role in the prediction of the 
branching ratios of neutral Higgs bosons. For example the decay mode 
h° — > bb could be significantly modified at large values of tan (3 and/or 
small values of sin a. The prediction of the branching ratios depends 
on the set of MSSM parameters, in particular the spectrum of SUSY 
particles change appreciably the SUSY Higgs decays. 



In Fig 4.2 we show the different decay modes of the neutral Higgs 
bosons for two different values of tan (3 as functions of the Higgs masses. 
The branching ratio of the charged Higgs boson is also shown in order 
to complete our analysis. As we know there is a limit for the light 
Higgs mass in the MSSM, m h o < 130 GeU[31], therefore h° will decay 
mainly into fermion pairs, in particular the most important channel is 
h° — > b b. This is in general also the dominant decay mode of the H° 
and A bosons, since for tan/5 ^> 1 the decay rate into bb and t~t + 
pairs are of the order of 90% and 10%, respectively. For large masses 
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Figure 4.2: Dominant MSSM Higgs bosons decay branching ratios as func- 
tions of the Higgs boson masses for tan/3 = 3 and 30. [from reference [34]] 



4.2. HIGGS BOSON DECAYS INTO TWO NEUTRALINOS 
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the top decay channels H°, A — > t t are suppressed for large tan (5. 

In order to complete our analysis, the SUSY decays of the neutral 
Higgs bosons must be considered, which could be dominant in differ- 
ent regions of the parameter space. In general any Higgs boson could 
decay into sfermions (squarks and sleptons), charginos or neutralinos. 
However taking into account the lastest results of the SUSY searches 
experiments [29], we know that in the case of the light Higgs boson, the 
only allowed SUSY channels are two neutralinos or two sneutrinos. For 
the heavy CP-even and CP-odd Higgs bosons the decays into squarks 
or charginos are also allowed, excluding the decay of A into two sneu- 
trinos. 

These SUSY decays will be dominant, of course, when the channels 
present in the Standard Model are suppressed. As we already noted 
the most important decay models are the decays into b or r pairs, these 
channels are suppressed in the case of low or moderate tan (3, or when 
the mixing angle a of the Higgs sector is quite small. Combining these 
two scenarios, we will able to get significant branching ratios for these 
channels. 

The various decay widths and branching ratios of the SM and MSSM 
can be calculated in a very precise way with the Fortran code HDECAY 
[39], where all the relevant experimental constraints are taken into ac- 
count. The subroutines of HDECAY dealing with the decays of neutral 
Higgs bosons into neutralinos use the results of reference [28]. 

4.2 Higgs boson decays into two Neutrali- 
nos 

In the Standard Model there are no invisible decays of the Higgs boson, 
since no v R is present in the model. In the MSSM the situation is quite 
different, there are new couplings which allow new decays of the Higgs 
bosons. The MSSM neutral Higgs bosons h°, H° and A could decay 
into invisible neutralino x? or sneutrino v pairs. These decay modes are 
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invisible in the case that the neutralinos or sneutrinos are the lightest 
SUSY particles (LSP) and -R-parity is conserved. Note that in the case 
of the CP-odd Higgs field, there is only one possibility the decays into 
two neutralinos, since the coupling to two sneutrinos does not exist. 
In this section we will describe in detail neutral Higgs decays into two 
neutralinos in the gaugino limit, considering quantum corrections at 
one-loop level. 

Before computing and discussing the partial widths for the decays of 
the neutral Higgs bosons into pairs of identical neutralinos, let us dis- 
cuss the properties of the couplings g% ~ ~ , where $ = h°, H° and A 

Xi Xi 

and X? is the lightest neutralino, in our case the lightest supersymmet- 
ric particle (LSP). 

At tree level, the couplings of the neutralinos x°i to the neutral CP- 
even Higgs bosons = h°, H° and to the CP-odd boson A° are given 
by: 

dlxX = 1 1^ 2 ~ tan^A^XtyV^ + e<p N j4 ) + i <-> j] (4.6) 

9%x°x° = I [(^2 - tan 9 w N a )(d A0 N j3 + e A oN j4 ) + i <-> j] (4.7) 
where the quantities d<$, and e$ are: 

d H o = — cos a , d h o = sin a: , d A o = sin/3 (4.8) 

e H o = sin a , e h o = cos a , e A o = — cos j3. (4.9) 

and Nij are the components of the matrix which diagonalizes the four 
dimensional neutralino mass matrix. 

Now using these equations we find the expressions for g% ~ - : 

#°°x?x? = 9 [( N i2 - t&n 6 w N n ) (sin aN 13 + cosaJV 14 )] (4.10) 
g° H o^o = g [(N 12 - tan e w N u )(- cos aN 13 + sinaN 14 )] (4.11) 
75 ^x?x? =^9 [(Ni2 - tan 9 W N U ) (sin [3 N 13 - cos/3JV 14 )] (4.12) 
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we see that all these couplings are exactly zero in the pure gaugino 
(JV13 = iV 14 = 0) or higgsino (N u = N 12 = 0) limit. 

Inserting eqs. (3.29) into eqs. (4.10) to (4.12), we see that the LSP 
couplings to the Higgs bosons $ = h°,H°,A° already receive contribu- 
tions at 0{m z ): 



9%xlx\ ~ d^N 13 +e^N u ~ s w m z 



fx 2 - Ml fi 2 - M 2 

(4.13) 

Similar expressions can be given for the couplings of the wino-like state. 

This suppression of the tree-level couplings follows from the fact that, 
in the neutralino sector, the Higgs boson couples only to one higgsino 
and one gaugino current eigenstate, together with the fact that mixing 
between current eigenstates is suppressed if ^> m z . These couplings 
thus vanish as \n\ — > oo. 

Knowing all the properties of the couplings, and taking into account 
that neutralinos are Majorana particles, we are able to compute the 
partial widths for the decays of the neutral Higgs bosons, $ = h°, H°, 
and A , into pairs of identical neutralinos: 

r(4>-xW) = ^fo ; + 9 i iK r. < 4 - 14 > 

where n = 3 for the CP-even fields h° and H°, while for the CP-odd 



4 m 2 



*0 , 



A Higgs boson n = 1. M$ is the Higgs mass, and /3$ 2 = (1 )■ 

We include the possible one-loop corrections to the couplings gi^o, 

Xl Al 

which will be discussed in the next section. 



4.3 Higgs decays into Neutralinos at one- 
loop 

At the one-loop level, the couplings of the lightest neutralinos to the 
Higgs bosons can be generated, in principle, by diagrams with the ex- 
change of either sfermions and fermions, or of charginos or neutralinos 
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together with gauge or Higgs bosons, in the loop. However, the latter 
class of diagrams can contribute to the couplings of Higgs bosons to 
neutralinos only if one of the particles in the loop is a higgsino. These 
loop contributions will thus be suppressed by inverse powers of 
in addition to the usual loop suppression factor, since the couplings 
g^H^XiXl) ~ 0(m w /fj), g (W ± x° 1 Xi ) ~ 0{m 2 w /^ 2 ) in the gaugino 
limit. We therefore do not expect them to be able to compete with the 
tree-level couplings that exist for finite 

We consider diagrams with fermions and sfermions in the loop, as shown 
below. For the couplings, only the third generation (s)particles, 

which have large Yukawa couplings, can give significant contributions 
to the amplitudes. Note that in the bino limit there is no wave func- 
tion renormalization to perform, since the tree-level couplings are zero. 
Off-diagonal wave function renormalization diagrams could convert one 
of the gaugino-like neutralinos into a higgsino-like state, but this kind 
of contribution is again suppressed by and can thus not compete 

with the tree-level coupling. 




The Feynman diagrams contributing to the one-loop couplings of the light- 
est neutralinos to the <3? = h°,H° and A Higgs bosons. Diagrams with 
crossed neutralino lines have to be added. 

We have calculated the contributions of these diagrams for arbitrary 
momentum square of the Higgs, finite masses for the internal fermions 
and sfermions as well as for the external LSP neutralinos, and taking 
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into account the full mixing in the sfermion sector. The amplitudes are 
ultra-violet finite as it should be. The contributions from diagrams a) 
and b) to the couplings are separately finite for each fermion 

species. 1 We have performed the calculation in the dimensional reduc- 
tion scheme [40, 41]; since the one-loop couplings are finite and do not 
require any renormalization, the result should be scheme independent. 
The results are given below for a general gaugino limit >> Ml, M 2 , 
for arbitrary ordering of Mi and M2). 

The one-loop Higgs boson couplings to the LSP neutralinos in the 
gaugino limit are given by: 



t>x\xl 



9 
4tt 2 



(/) 



(4.15) 
(4.16) 



where 



{s2e f (vi + 



v 3j 



m) +m?o)C (/i)+4m?oCi + (/i 



m f9<t>ff 

+ (m\ +m) + m| ; ) C (/ 2 ) - 4m|oC 1 + (/ 2 )^ 
+ 2{ Vl + v 2 c 2 d} )m f m^ [c (/i) - 2Ci + (/i/ 

+ 2( Vl + v 2 s 2 ef )m f m^ [C (/ 2 ) - 2C+(/ 2 )] } 

C^/i {-^-(^1 + v 3 )m f C (fi, fi) + 2(^i + v 2 c 2 e _)m ii0 Cf(f 1 , A) J 
^0/2/2 { s 29 7 (ui + v 3 )m f C (f 2 , h) + 2(^i + v 2 s 2 e} )m^Cf{f 2 , / 2 ) j 
% i/2 {-2c 20/ >i + U3 )m / C (/i, / 2 ) - 2s 2 ^ 2 m^C'i + (/i, / 2 )} (4.17) 



_ m f9A0ff f ( v 



(m^ - mj - m|o) C (/i) - (m| - raj - ra|o) C (/ 2 )] 



1 The contribution of diagram a) is finite only after summation over both sfermion 
mass eigenstates. 
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m f g A o ff 



|2K + v 2 c 2 e .)m^m f C (f 1 ) + 2( Vl + v 2 s 2 ef )m^m f C (f 2 ) 



+ C A°hh {~ 2m f( v i + V 3)C U\, f 2 ) + 2m^>(2v 1 + v 2 )s 2df C 1 (/i, / 2 )} 

(4.18) 

The Higgs-fermion-fermion coupling constants are given by 



9h°uu — 

9h°uu — 



cos a 
sin P 
sin a 



sin (3 
9a°uu = cot f3 



9h°dd — 

9H0dd = 



sin a 

cos (3 
cos a 



cos/3 ' 
9A»dd = tan/? , 



(4.19) 



and 






Vi = 


\ (gQ f N n tan e W) 2 , 


v 2 = 


2Q/ 


-2 + -f + 2 


v 3 = 







J 3 ' \ jV- 



12 . lL(. 



N- 



12 



Qf J N ±1 tan 9 W Qf \N U tan6 w 



1 . 



(4.20) 



while the Higgs-sfermion-sfermion coupling constants read: 



™ 2 /-i 2 

J 3 c a ~ ^uS w C 2 g- 



m 2 u g h o uu m u s 29 . 



mw 2m 



w 



[A u gh°v.u + H9h° 



n - mz 

^h°U2U2 — s /3+c 



^3 S 9 a + QuS w C 2 g. 



m u g h o uu m u s 26 . 

+ — y%i9h°uu + H9h°i 



mw 2m 



w 



<^h uiu 2 — s /3+c 
CW 



Q u s w - 1^/2 s 2e - - [A u g h o uu + ng H o U u\ c 26> fi 

J zmw 
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n - mz 

U h°d idi - — s /3+c 

cw 



T t c l d _ QdSwC 2 d 



mjg h o dd m d s 2 e d 
mw 2mw 



[Adgh°dd - H9H°dd\ 



m z 



cw 



Ch°d 2 d 2 - — s i3+a I^s 2 ,. + QdS w C 2 e 3 



v./.- . "' dS -"> [A d g h o dd - ng H o dd ] 



m w 2m w 



^h"d!d 2 ~ S l 3 + a \Qd s W ~ 

Cw L 



m d 



$2e s - [A d g h o dd - fj>g H »dd] C26 d 



Cj? Ml«l 



m z 
cw 



tu 2 /"i 'z 
1 3 % ~ { 4uS w C 2 8 il 



m u g H o uu m u s 2e ~ 



mw 2m 



w 



[A u gH° uu ~ ^9h°i 



Ch°u 2 U2 ~ 



m z 
cw 



I3 S L + QuSwCWz 



m 2 u g H o uu m u s 2e . 

H — ~ [ a u9h°uu — H9h 



mw 2m 



w 



Cfl»uii 2 — -cg +a [QmS^ — -^3/2 S2e c — - — — [A u gH°uu — l^9h°uu\ c 2e n 

cw zmw 



m„ 



C 



H°didi 



m z 
c w 



-C/3+a 



lie 2 * - Q d s w c 2e , 



m d g H °dd m dS 2 e A 



mw 2m 



w 



[A d g H °dd + fJ>9h dd} 



C 



H^d 2 d 2 



m z 
c w 



C/3+a 



+ QdS w c 2Q - 



m z d g H o dd m d s 26 - 

+— [A d gH°dd + l^gh°dd\ 



mw 2m 



w 



r - mZ n 

u H0d!d 2 — „ c P+a 
Cw 



m d 



QdS w - J3/2J s 2 e d - - [AjffflOdd + l*9h°dd] C28 d 



m„ 



Caos iS2 = — - — (A u cot/3 + /x) 



2m 



m d 



(4.21) 
(4.22) 
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with s a+ p = sin(o; + f3), c a+ p = cos(a + (3), s$. = sin 6^, eg. = cos9j 

and = (— )|. The Passarino-Veltman three-point functions, are 
defined as 

Cof(f) = C +r(g 2 ,m|o,mJ,m5,mJ-); 
CSi + '-(/i,/ 2 ) = C^-(q 2 ,ml 0i ,ml,ml,m}), (4.23) 
see Appendix A for the explicit form of these functions. 

Now knowing all the details about the tree level couplings and the 
quantum corrections considered, we are able to show few numerical ex- 
amples to see the effect of the quantum corrections to the branching 
ratios of the neutral Higgs bosons decays into two neutralinos. 

The real part of the coupling of an on-shell lightest h° boson to a LSP 
pair is displayed, at the Born and one-loop level, in the Figs. 4.3 and 
4.4. It is shown as a function of the sfermion masses, for tan/5 = 15, 
\x = 1 TeV and pseudoscalar mass input values M A o = 200 GeV 
(Fig. 4.3) and 1 TeV (Fig. 4.4). Top quarks couple with 0(1) Yukawa 
coupling to the Higgs bosons, and for the given choice of large \A t \ the 
(dimensionful) h°titi coupling significantly exceeds the t\ mass. More- 
over, due to tL — tjt mixing the lighter t mass eigenstate is often not only 
lighter than the other squarks, but also lighter than the sleptons. If \A t \ 
is large, as in the present example, the loop corrections to the h° Xi Xi 
coupling can even exceed the tree-level contribution. The variation 
of the one-loop contribution to the coupling is again mostly due to the 
natural decrease with increasing masses of the sfermions running in the 
loop, which decouple when they are much heavier than the h° boson. 
For rriq ~ 420 GeV [i.e. raj ~ 210 GeV], is near its experimental 
lower bound of ~ 100 GeV, due to strong II —Ir mixing. This implies 
that m~ tx will grow faster than linearly with increasing rriq, which ex- 
plains the very rapid decrease of the loop corrections. However, there 
is also a variation of the tree-level coupling for M40 = 200 GeV which, 
at first sight, is astonishing. It is caused by the variation of the mixing 
angle a in the CP-even Higgs sector, and to a lesser extent by the vari- 
ation of M^o, due to the strong dependence of crucial loop corrections 
in the CP-even Higgs sector on the stop masses. In fact, for the set 
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Figure 4.3: The lightest h° boson couplings (top) and branching ratios 
(bottom) to pairs of the lightest neutralinos as functions of the common 
slepton mass. These results are given for tan f3 = 15, fj, = 1 TeV, rrig = 
2mj, A t = 2.9mg and gaugino masses M\ = 30 GeV, M2 = 120 GeV, and we 
took M40 = 200 GeV 
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Figure 4.4: The lightest h° boson couplings (top) and branching ratios 
(bottom) to pairs of the lightest neutralinos as functions of the common 
slepton mass. Most parameters are as in Fig. 4.3, and we took M A o = 1 TeV 
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of input parameters with M^o = 200 GeV at small slepton masses, we 
are in the regime where sin a, which appears in the h°XiXi coupling 
[and which enters the h°bb coupling as will be discussed later], varies 
very quickly. This "pathological" region, where the phenomenology of 
the MSSM Higgs bosons is drastically affected, has been discussed in 
several places in the literature [42, 43]. 

The branching ratios for the decays of the lightest h° boson are shown 
in the Figs. 4.3 and 4.4, for the same choice of parameters previ- 
ously discussed. They have been calculated by implementing the one- 
loop Higgs couplings to neutralinos in the Fortran code HDECAY [39] 
which calculates all possible decays of the MSSM Higgs bosons and 
where all important corrections in the Higgs sector, in both the spec- 
trum and the various decay widths, are included. The branching ratio 
BR(/i° — > XiXi) can aheady exceed the one permille level with tree- 
level couplings. After including the one-loop corrections, the branching 
fraction BR(h° — > XiXi) can be enhanced to reach the level of a few 
percent. 

The branching ratio is especially enhanced if the usually dominant de- 
cay into bb pairs is suppressed, i.e. if | sino;| is very small; recall that 
the h°bb coupling is oc sin a/ cos (3. In our examples this happens for 
M^o = 200 GeV and m- { ~ 250 GeV. In this case fi^o^o is only about 
half as large as in the decoupling limit M40 — > oo, but the loop contribu- 
tion to this coupling is still sizable for this value of the sfermion masses, 
and has the same sign as the tree-level coupling, leading to a quite large 
total coupling. The branching ratio falls off quickly for smaller sfermion 
masses, since here sin a becomes sizable (and positive). Moreover, for 
mi ~ 210 GeV the tree-level and one-loop contributions to the cou- 
plings have opposite sign. The branching ratio also decreases when mi 
is raised above 250 GeV, albeit somewhat more slowly; here the rapid 
decrease of the loop contribution is compensated by the increase of the 
tree-level coupling, which however does not suffice to compensate the 
simultaneous increase of sin 2 a. 

Fig. 4.5 shows the dependence of the h° Xi Xi coupling and of the 



46 CHAPTER 4. SUSY DECAYS OF NEUTRAL HIGGS BOSONS 



corresponding h° branching ratio on the mass of the LSP, m^o ~ Mi. 
We see that the tree-level contribution to this coupling depends essen- 
tially linearly on the LSP mass. Eqs. (4.10) and (3.29) show that, for 
the given scenario where c Q ~ sp ~ 1, this linear dependence on Mi 
originates from JV14, where the contribution with Mi in the numerator 
is enhanced by a factor of tan /3 relative to the contribution with /i in 
the numerator. Therefore the contribution oc Mi is not negligible even 
though in Fig. 4.5 we have Mi <C \fi\. On the other hand, the one-loop 
contribution to this coupling depends only very weakly on Mi. The 
small increase of this contribution shown in Fig. 4.5 is mostly due to 
the explicit m^o dependence of the loop coupling (equation 4.17); the 
change of N u with increasing Mi, as described by eq. (3.29), plays a 
less important role. The increase of the total coupling with increas- 
ing Mi nevertheless remains significant. However, Fig. 4.5 shows that 
for m^o > 15 GeV this increase of the coupling is over-compensated 
by the decrease of the f3 z threshold factor in the expression for the 
— >■ XiXi) partial width. 

Once one-loop corrections are included, for certain values of the MSSM 
parameters the branching ratio for invisible h° boson decays can thus 
reach the level of several percent even if Xi is an almost purely bino. 
This would make the detection of these decays possible at the next gen- 
eration of e + e~ linear colliders. At such a collider it will be possible to 
isolate e + e~ — > Z°h° production followed by Z° — > l + £~ decays (£ — e 
or fi) independent of the h° decay mode, simply by studying the distri- 
bution of the mass recoiling against the £ + £~ pair. This allows accurate 
measurements of the various h° decay branching ratios, including the 
one for invisible decays, with an error that is essentially determined by 
the available statistics [44]. Since a collider operating at y/s ~ 300 to 
500 GeV should produce ~ 10 5 Z°h° pairs per year if | sin(o; — f3)\ ~ 1 
one should be able to measure an invisible branching ratio of about 3% 
with a relative statistical uncertainty of about 2%. 

We now turn to the heavier MSSM Higgs bosons H° and A . The 
couplings to the lightest neutralinos are shown in Figs. 4.6 and 4.7, for 
the same input parameters as in Fig. 4.3. As can be seen, up to a rela- 



4.3. HIGGS DECAYS INTO NEUTRALINOS AT ONE-LOOP 47 









_tg|3=15/M,=200GeV 






( u=1TeV/M.,. pl =250GeV 






— one — loop level 






— tree level y 













1 10 20 30 40 50 

M, [GeV] 







_ BR(h^ x ; x ?) \ 




. tg^=15/M A =200GeV \ 




. M=1TeV/M,„ pl =250GeV \ 




_ — one-loop level \ 




- ■ - tree eve 









1 10 20 30 40 50 

M, [GeV] 



Figure 4.5: The lightest h° boson couplings (top) and branching ratios (bot- 
tom) to pairs of the lightest neutralinos as functions of Mi . The parameters 
are as in Fig. 4.3 with mj = 250 GeV and M A o = 200 GeV. 
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Figure 4.6: The heavier CP-even Higgs boson coupling H (top) and 
branching ratios (bottom) to pairs of the lightest neutralinos as functions of 
the common slepton mass. The parameters are as in Fig. 4.3. 
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Figure 4.7: The pseudoscalar A boson couplings (top) and branching ratios 
(bottom) to pairs of the lightest neutralinos as functions of the common 
slepton mass. The parameters are as in Fig. 4.3. 
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tive minus sign, the tree-level couplings of these two Higgs bosons are 
approximately the same since we are in the decoupling regime where 
cIao — —d H o, and e^o ~ e#o with |e^o| <C 1, see eq. (4.8). Eqs. (4.10- 
4.8) and (3.29) also show that the tree-level couplings of the heavy 
Higgs bosons exceed that of the light Higgs boson h° by a factor tan/3/2 
[ignoring contributions to eqs. (3.29) with Mi in the numerator]. On 
the other hand, the loop corrections are smaller in case of the heavy 
Higgs bosons. The corrections to the H° Xi Xi coupling are reduced 
by about a factor of 2 compared to the corrections to the h° xi xi 
coupling, mostly due to the relatively smaller coupling to t pairs, see 
eqs. (4.3). 

The corrections to the A xi Xi coupling are even smaller, since the 
CP-odd Higgs boson A cannot couple to two identical squarks. The 
contribution with two t\ squarks and one top quark in the loop, which 
dominates the corrections to the couplings of the CP-even Higgs bosons 
for small mj, does therefore not exist in case of the A boson. As a re- 
sult, the corrections to the coupling of A are not only smaller, but also 
depend less strongly on mf recall that for our choice of parameters 
increases very quickly as mj is increased from its lowest allowed value of 
~ 210 GeV, which comes from the requirement > 100 GeV. Note 
also that the H°tt and A°tt couplings are suppressed by a factor cot f3 
relative to the h°tt coupling, see eqs. (4.19); this becomes important 
for large squark masses, where the one-loop corrections are relatively 
less important. Altogether we thus see that the one-loop corrections 
are much less important for the heavy Higgs bosons. Note also that 
for the given set of parameters they tend to reduce the absolute size of 
these couplings. 

Again, because in the decoupling regime the CP-even H° boson and 
the pseudoscalar A boson have almost the same couplings to Standard 
Model particles and to the neutralinos [at the tree level], their branch- 
ing ratios are approximately the same. The one-loop contributions 
decrease the branching ratios by at most ~ 10 to 40%. Note that the 
total decay widths of the A and H° bosons are strongly enhanced by 
tan 2 /3 factors [T(H°, A — > bb) oc m%g\ #o & &]■ This over-compensates 
the increase of their couplings to neutralinos, so that their branching 
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ratios into Xi pairs are far smaller than that of the light Higgs boson h°, 
remaining below the 1 permille level over the entire parameter range 
shown. Moreover, the cross section for the production of heavy Higgs 
bosons at e + e~ colliders is dominated by associated H°A° production, 
which has a much less clean signature than Z°h° production does. 
Branching ratios of the size shown in Fig. 4.6 and 4.7 will therefore not 
be measurable at e + e~ colliders. In fact, they will probably even be 
difficult to measure at a collider "Higgs factory" ; recall that the 

Z° factories LEP and SLC "only" determined the invisible decay width 
of the Z° boson to ~ 0.1%. 

After our analysis was completed a different group studied the one-loop 
corrections to neutral Higgs bosons decays into neutralinos in the gen- 
eral case, where we have the decays Hf — > XmXn (i=l,2,3) [45]. They 
confirmed our results computing the branching ratios for similar values 
of the parameters. 

As we mentioned before we consider all our parameters to be real, 
however in the general case there are new CP violating phases in the 
MSSM. Due to the potentially important role of the decays analyzed 
above, reference [46] studied the dependence of the branching ratios 
on the CP violating phases. Moreover a large correlation between the 
spins of the x states produced in the decays of heavy neutral Higgs 
bosons was found. 



52 CHAPTER 4. SUSY DECAYS OF NEUTRAL HIGGS BOSONS 



Chapter 5 



The Minimal 
Supersymmetric SU(5) 

5.1 SUSY and Unification of the gauge 
couplings 

In the Standard Model for each group a gauge coupling constant is 
defined. It has been known for a long time how the gauge couplings 
change with energy [47]. If we study the evolution of the gauge cou- 
plings we see that in the context of the Standard Model the gauge 
couplings never meet. 

The meeting of the gauge couplings in the Minimal Supersymmetric 
Standard Model is an impressive prediction [2, 3, 4, 5] (see Fig. 5.1), 
which tells us that at the high scale M GUT ~ 10 16 GeV, all the inter- 
actions are unified. Above the GUT scale the gauge couplings remain 
together only if new particles are present, this is the case of SUSY 
SU(5). 

5.2 Particle assignment 

The minimal Supersymmetric SU (5) model is the simplest framework 
where the unification of the Standard Model Interactions is realized. 



53 



54 



CHAPTER 5. THE MINIMAL SUPERSYMMETRIC SU(5) 




Figure 5.1: Evolution of the inverse of the three coupling constants in 
the Standard Model (left) and in the super symmetric extension of the SM 
(MSSM) (right). Only in the latter case unification is obtained. The SUSY 
particles are assumed to contribute only above the effective SUSY scale 
A-SUSY of about 1 TeV, which causes a change in the slope in the evolution 
of couplings. The thickness of the lines represents the error in the coupling 
constants as of 1991 [48]. 



5.2. PARTICLE ASSIGNMENT 



55 



Using the quantum numbers of the SM particles Georgi and Glashow 
[49] showed how the matter is unified partially in two irreducible repre- 
sentations 5 and 10. Using the SU(3) x SU(2) decomposition of these 
representations, the fermions of one family are accommodated as: 
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SU (5) is a Lie group of rank 4, with 24 generators. Therefore we will 
have 24 gauge fields in our model, the usual Standard Model gauge 
bosons plus 12 additional gauge bosons: 



A M (24) = -A°A», a = 1,2... 24. 



(5.1) 
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where the A a are given by: 
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The SU(3) x 577(2) decomposition of the 24-plet is given by: 
24 = (8, 1)0(3, 2) 0(3, 2) 0(1, 3) 0(1,1) 

A(24) = G y 0(X,, Fj) ©fa X,) 0(Vy + , VU 3 , W") B c 



The £77 (3) octet Gij are identified with the gluons. The SU(2) doublet 
(Xi,Yi) represents the two superheavy SU(3) triplets X and Y gauge 
bosons with electric charges 4/3 and 1/3 respectively. The SU (2) triplet 
(W + , W 3 , W~) are identified with the SM 577(2) vector bosons and fi- 
nally B° is the U(l) vector boson. 

Now in order to know if our model reproduces the well known low- 
energy physics, we must break the symmetry spontaneously to the 
Standard Model gauge group SU(3) C x SU(2) L x U(1) Y - In order 
to achieve this we define the minimal Higgs sector, which is composed 
of three representations, 5h, 5h and the adjoint representation £(24): 



T 3 

Hi 

H° 2 



\ 



Ti 

v -m j 



S(24) 



-(3,2) 



-(3,2) 




J 24 



Knowing all the particles of our model we are ready to write down 
the interactions and analyze the possible predictions coming from new 
interactions. 



5.3 The SU(h) Lagrangian 

Using the tools given in Chapter 1 and introducing a superfield for 
each representation of SU (5), we can write the lagrangian of our model 
[50, 51]: 
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4S = i J d 2 9 Tr(W«W 5 , a ) + - J d 2 9 Tr(W«W 5 , a ) 
+ ]T J d 2 9 d 2 9 e ^^(24) ^ 

$=10,S ^ 

+ jd 2 9W 5 + j d 2 9 W 5 (5.2) 

The most general (in the renormalizable limit) superpotential W5 of 
SU(5) which is R-parity invariant, W5 has two important pieces, the 
corresponding to the Higgs self-interactions and other describing Yukawa 
couplings: 

W 5 = W H + Wy (5.3) 
The superpotential of the Higgs sector reads as: 

W H = TrS 2 + ^ TrS 3 + 77 5*115* + 771*5*5* (5.4) 
while the Yukawa superpotential is: 

W Y = 10IV105* + 10T D 55 H (5.5) 
where T's are 3x3 Yukawa matrices. 

In the supersymmetric standard model language the Yukawa sector 
can be rewritten as 

W Y = QY V U C H + HQY D D C + HLY E E C 

+ QAQT + U C BE C T + QCLT + U C DD C T (5.6) 

where except for the heavy triplets T and T the rest are the MSSM su- 
perfields in the usual notation. The generation matrices Y\jp,E an d A, 
B_, C_ and D_ can in general be arbitrary. In the minimal SU(5) defined 
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above one finds A — B_ — Y v — Yy = Tjj, and C_ = D_ = Yd = Yj = 
at the GUT scale. 

From these interactions we can find the different effective operators 
contributing to the decay of the proton. These are LLLL and RRRR 
operators: 



^ J d 2 9 (Q AQ)(QC L) (5.7) 
J d 2 9 (U c B E c ) (U c D D c ) (5.8) 



Mi 



In the next Chapter we will study all the properties of these operators, 
and we will analyze the predictions in the minimal model. 

5.4 Symmetry Breaking 

We need the following symmetry breaking: 

SU(5) x SUSY =^ SU(3) C x SU(2) L x U(1) Y x SUSY 

To study this we have to use W# and calculate the relevant F-terms and 
set them to zero to maintain supersymmetry down to the electroweak 
scale. Computing the F-terms and using the condition TrS = we find 
the possible solutions which give us the symmetry breaking preserving 
SUSY: 

Case 1. 



(E) = (5.9) 
In this case the SU (5) symmetry remains unbroken. 

Case 2. 
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In this case S77(5) breaks down to S77(4) x U(l), and the last possible 
solution is: 



Case 3. 



< E > = x 
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(5.11) 



This is the desired vacuum since 577(5) is broken to SU (3)c x SU (2)l x 
U(1)y- In the supersymmetric limit all vacua are degenerate. To com- 
plete the symmetry breaking, the Gsm must be broken to SU(3)c x 
U(l) em . This is caused by the following expectation values: 



<5h> 



/ \ 




\ V2 J 



(5h) = 



V 



o 





-JiL , 
V2 / 



(5.12) 



, is 
Y T 



The fact that from W Y we get A = B = Y v , C = Y E , D = Y D 
simply a statement of SU(5) symmetry. On the other hand Yy = 
and Yd = Y^ result from the SU(4) Pati-Salam like symmetry left 
unbroken by (5#) and (5#). Under this symmetry d c <-> e, u <-> u c , 
d <-> e c . Of course, this symmetry is broken by 7^ where 
a = 1,2,3; this becomes relevant when we include higher dimensional 
operators suppressed by {T,)/M Pl , which will be considered in the next 
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section. 

Knowing how SU(5) is broken to the Standard Model gauge group, 
we can compute the masses of different Higgs superfields in our theory. 
From the expression of Wh we can compute the Triplet and the \i mass: 

M T = 2rj^ + m H (5.13) 
A 

and 

^ = m H -3r]^y (5.14) 

The triplet mass Mt must be close to the GUT scale, while \i must 
be close to mw- From the above relations we see that only when the 
parameters of the potential are fine-tuned, we can explain this differ- 
ence in the masses. It is the Fine Tuning problem of GUTs. As we 
see here we need a lot of fine-tuning to explain how the Triplet is much 
heavier than the doublet, this is the so-called Doublet-Triplet Splitting 
or Hierarchy Problem. Supersymmetry only helps us to stabilize the 
splitting againts radiative corrections, but by itself does not explain its 
origin. 

When the symmetry is broken the X and Y gauge bosons become 
massive, 

Mx = M Y = 5^g 5 ^y (5.15) 
while we find for the members of S 

5 

m 8 = m 3 = -mj (5.16) 
2 

Note that in this case the weak triplet and color octet masses are equal. 



5.5 Fermion masses 



As we mentioned before in the Minimal Supersymmetric SU (5) model 
we find the relation Y D = Yj = T D at the GUT scale. When 5# gets 
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the expectation value (5#)= diag ((J, 0, 0, 0, — the quark and lepton 
masses are related as: 



m e = m d , m^ = m s , m T = m b (5.17) 

The first two relations are wrong. Note that the incorrect relation 
— = ^ is predicted to be valid at any scale. On the other hand the 
relation for the third generation can be considered a great success of 
the theory. 



We can imagine many ways to improve the mass relations for the first 
two generations [51], but the simplest and most suggestive one is to 
include l/Mpi suppressed operators which are likely to be present; af- 
ter all due to the small size of the Yukawa couplings, these operators 
should be more important for the first two generations where the theory 
fails, and they require no change in the structure of the theory. Note 
that the value of the ratio M $£f ~ 10~ 3 — 10" 2 is even bigger than the 
Yukawa couplings of the first generation. 

The explicit form of the renormalizable, and all the relevant non-renormalizable 
terms are 



(ym ym 
lWUlflSS + l^f la M l ^ n H + 10^10^-/- 
Mpi IVlpi 

+ 5 Hi 10^g ab 5 bj + 5 Hl ^W a k g lab 5 bk + 5 Hi 10% g 2ab ^5 bk (5.18) 

Mpi Mpi 

where i, j, k,l,m,n are SU(5) indices, and a,b — 1,2,3 are generation 
indices. 

After taking the SU(5) vev (S) = a diag(2, 2, 2, -3, -3) we get at 
M G ut scale. 
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M P ,y "> m pi 

£ = - 9 - 2 £ 91 - 2 ^ < 5 - 19) 

using A a — m-s. 

Note the relation Y E -Y D = C- D 

In the limit Mpi — > oo we recover the old relations, but for finite 
a/Mpi ks 1CT 3 — 1CT 2 one can correct the relations between Yukawas 
and at the same time have some freedom for the couplings to the heavy 
triplets. 

Clearly, due to SU(5) breaking through (E), the T, T couplings are 
different from the H, H couplings. However, under the SU(4) symme- 
try discussed before A <-> B_, C <-> D and Y v <-> Fj\ Only the terms 
that probe (E£) — (E|) can spoil that; this is why /x and g\ still keep 
Yu = YT,A = B and C — D. 

Now we are ready to improve the mass relations for the first two fam- 
ilies. In order to get the correct relation ^ ~ \—, we must impose 



specific values to the coupling g 2 since: 



If we assume that Yd and Ye are diagonals, and using the relations 
m e = ^ and m M = 3m s , g 2 must satisfy the following relations: 



5.6. sm 2 e w 
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(ftOn = -~ m A Mpi (5-21) 
vy2ni 15 (if) m s V ; 

and 

Note how the parameters of the scalar potential enter in the expressions 
for fermion masses. 



5.6 sin 2 6^ 

In the Standard Model the Weinberg angle 9w is a free parameter, 
which plays an important role in weak interactions. From experiment 
we know that sm 2 9 wl ^(M z ) = 0.23117 ± 0.0016 [53]. 

Let us see what happens in SU(5). In the Standard model we define 
tan#yK = 9u(1) , the electromagnetic charge operator Q em = T 3W + ^ 

must be part of the £77(5) operators, therefore TrQ em = 0. 

Now using for example the fundamental representation 5h we can pre- 
dict QemiT) = — |, at the same time using 5 representation we can 
predict Q em (d c ) = ^Q em (e). This is one of the most beautiful predic- 
tions of GUTs, the quantization of the electric charge. 

For any fundamental representation we get: 

1 \ 

10 

10 

-3 

V J 



Qem(5) — — - 



66 



CHAPTER 5. THE MINIMAL SUPERSYMMETRIC SU(5) 



A23 
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-1 J 



knowing these two operators, we see that the hypercharge operator 
must be: 



-(») 

2 K 1 
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as in SU (5) all the couplings are equal, we can get the relation between 
gu(i) an d 9su{2)- From the relations listed above we can write the 
following expressions: 



Y 



\24 



(5.23) 



T: 



9su(2)- 



3W 



9b- 



\23 



(5.24) 



or 



concluding that ^(i) = yjfg 5 and g su{2 ) = 95, so tan 2 9 W = 
sin 2 6w — §■ It is one of the most important predictions of super- 
symmetric Grand Unified Theories and in particular of SUSY SU(5). 
Note that this value is at the GUT scale, when we use the renormal- 
ization group equations and compute the value of this quantity at the 
electroweak scale, we see that it agrees with the experimental measure- 
ments [54, 55]. 



Chapter 6 

Proton Decay in the 
Superworld 

6.1 B violating operators 

As we know the Baryon (B) and Lepton (L) numbers are conserved 
in the Standard Model. It is a consequence of the particle assignment 
and the gauge principle. However these symmetries could be broken 
at a high scale M ^> mw- If at the high scale this happens, then we 
will have an effective operator C e ff which describes the new possible 
interactions: 

Qd 

where O d represents an operator of mass dimension d, Cd is a coefficient, 
and D is the space-time dimension. Note that in our "rea/" world we 
have D = 4. 

If the Baryon number is broken, we have a new prediction, the decay 
of the proton. This is the case of grand unified models such as SU(5) 
and 5*0(10), where from the matter unification we get new effective 
operators of the type 6.1 [56, 57]. In the case of non-conservation of 
Leptonic number, we have the possibility to explain the smallness of 
neutrino masses, due to the presence of a large Majorana mass term 
[58, 59, 60, 61]. 
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Using the superfields of the Minimal Supersymmetric Standard Model 
[Q = (U,D), L = (N, E), U c , D c , E% we can write down all 
the possible effective operators contributing to the decay of the proton, 
which are SU(3) C x SU(2) L x U(1) Y invariant [62, 63, 64, 65, 66, 67]. 

d = 4 operators 1 : 




(6.2) 




(6.3) 




(6.4) 



d = 5 operators: 




(6.5) 




djkl 



(6.6) 



d = 6 operators: 




(6.7) 




(6.8) 



Note that these are the operators present in Wnr (see equation 3.21) 
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where a, (3 and 7 are color indices; m, n, p and q isospin indices, while 
i, j, k and 1 represent generation indices. 

Note that the product of two d = 4 operators and d = 6 operators 
lead to proton decay at tree level. In the case of the dimension d — 4 
operators we have contributions with two fermions and one scalar field, 
the exchange of the scalar field can mediate proton decay. The d = 6 
operators directly yield terms with four fermions contributing to the 
decay of the proton. However d — 5 operators are quite special, in 
each term we have two fermions and two scalars fields, therefore they 
contribute only at one-loop once we dress these operators. 

In the case that extra spacetime dimensions are considered we see that 
there will be many new contributions to proton decay, without the 
suppression factor [68]. This is in our opinion the most important phe- 
nomenological problem of many models with extra dimensions. 

Knowing all the possible operators contributing to the nucleon decay, 
the general expression for the proton lifetime could be written as: 

r p oc \A d=4 + A d=5 + A d=e \- 2 (6.9) 

where Ai are the different amplitudes. Note that in four dimensions A4 
does not have any suppression factor, therefore the coefficients related 
with these contributions must be very small or maybe it is more natural 
find a symmetry to forbid these operators. 

The second possibility is realized, if we introduce a symmetry called 
Matter Parity (see section 3.6): 

M = (-l) 3 ^) (6.10) 

where M — —1 for Q, L, U c , D c , E c and M = 1 for H, H and G M . 
If we assume that this symmetry is conserved, we remove the d — 4 
contributions, retaining the d — 5 contributions as the most important 
ones. 
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Note that relation between R and M parities, R = (— 1) 2S M. The 
case where the R parity is not an exact symmetry has been analyzed 
in reference [69]. However the most interesting case is when i?-parity 
is conserved. As we mentioned before in this case we have an ideal 
candidate to describe the Non-Baryonic Dark Matter present in the 
Universe. Also the conservation of this symmetry is predicts in a large 
class of Grand Unified Theories as Minimal SUSY SO (10) [32]. 

6.2 d = 6 operators 

Let us analyze in detail the d = 6 contributions. We can use as an 
example the hermitian of the operator 6.7, working in 4 dimensions, 
choosing m — 1, n — 2, and using the properties of the Grassmannian 
variables we find the following four fermions effective operator: 



therefore we will have a contribution to proton decay at tree level, in 
this case the proton decay into tc° and e + , the usual most important 
channel coming from the d=6 contributions in grand unified models. 



1 




(6.11) 



u 



u 




d = 6 contribution to the decay of the proton. 



Usually the d = 6 processes are mediated by new superheavy gauge and 
Higgs bosons present in grand unified models. There are many aspects 
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to be considered when we compute the proton decay amplitudes. In 
the first place these operators are given at the GUT scale, therefore to 
compute the values of the lifetime, we must compute the matrix ele- 
ments for each channel, and study the evolution of these operators to 
the proton mass scale 1 GeV (see [70, 71, 72] for more details). 

Also there is a very important point related with the fermion masses 
and the prediction of proton decay. Assume that the fermion mass 
matrices are diagonalized as: 

U T Y U U C = Y* , D T Y D D C = Y* , E%Y%E = Y* , (6.12) 
when u — > U u, u c — * Uc u c and so on. 

Now if we write the d — 6 operators (eq. 6.11) in the physical ba- 
sis, we get: 

^ (e c y (E c D) d (u c y (U c U) u (6.13) 

As we see in general these operators depend of the textures for Yu, Yd 
and Ye, this means that the proton decay predictions will be different 
in each model for fermion masses [73, 74]. 

In the Minimal Supersymmetric SU(5), where from the relation of the 
mass matrices Yjy = Yy and Yd = Y^ we have Uc = U and D = Ec, 
the d = 6 operators are independent of the explicit form of textures, 
however as we mentioned in the last chapter this is not a realistic case, 
due to the problem of the mass relation for the first two families. 

We could say that d = 6 proton decay provides a way to test models of 
fermion masses, however as we know Mqut ~ 10 16 GeV which gives us 
r p (d = 6) ~ 10 35 years, a value which is much bigger than the present 
experimental bounds [75]. Therefore it is difficult to test the d = 6 
predictions at present experiments. See [76, 77] for the predictions in 
string-derived models. 
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6.3 d = 5 operators 

The d — 5 contributions are the dominant to the proton decay. They 
are mediated by the superpartner of the colored Higgses (Triplets), 
which are present in SUSY GUT models. In this case we have only -jg- 
as suppression, where is the Triplet mass. 

Let us understand how from these operators we get the proton decay 
amplitudes. We can use the following operator as a example: 



^ T.T.I 



LLLL 



M 7 



J d 9 € rnn dpq 6 a /3y Qima Qjn/3 Qkp"f L\q (6.14) 



in this case using the properties of the Grassmannian variables we see 
that for each contribution, there are two fermionic and two bosonic (or 
superpartners) fields. For example if % — j — I — 1 and k = 2, we find 
the following contribution to the decay into K + and u: 



-^fLLLL Udsu (6.15) 

Now we must dress this operator to find the four fermions operator 
contributing to proton decay. This is possible using gauginos and hig- 
gsinos, since these are Majorana particles: 



u u 




d = 5 contribution to the decay of the proton. 

From this graph we can appreciate that these contributions are present 
at one-loop level, where we have superpartners inside the loops, the 
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loop factor and the suppression must be considered to compute 

the proton decay amplitude. As we mentioned in the last section, these 
operators are valid at the GUT scale, therefore we must compute the 
matrix elements and study the running down to 1 GeV. However in this 
case there are many new factors to be considered. is usually the 
Triplet mass, which could be smaller than the GUT scale, therefore we 
must compute this in order to estimate the amplitudes. 

In the case of d = 6 operators we showed how the amplitudes could 
depend on the textures of fermion masses. For the d — 5 contributions 
we see that there is something new, the sfermion masses also appear in 
this case. Therefore in general the d — 5 contributions will depend on 
the textures for sfermions and fermions, since in a general SUSY model 
these textures are different. 

As we see in order to compute the proton decay amplitudes we must 
consider many unknown factors: the loop factor, which depends on 
the SUSY spectrum and mixings between fermion and sfermions, the 
Triplet mass and the matrix elements. Therefore we can conclude that 
it is very difficult to test the SUSY GUT models using proton decay, 
since the dominant d — 5 contributions are quite model dependent. 



6.4 Proton decay in Minimal SUSY SU(5) 



In the last chapter we studied the structure of the Minimal Supersym- 
metric SU (5) model. We noted that new interactions which violate the 
baryon and lepton numbers are present, when the matter unification 
is realized. From these new interactions we find the d — 5 and d = 6 
operators contributing to the decay of the proton. 

As we mentioned above the most important contributions are those 
with d — 5. From the superpotential of £77(5) we find the LLLL and 
RRRR d — 5 operators, which read as: 




(6.16) 
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±-B^D kl J <?6{U? Ef) (US Df) (6.17) 

Knowing these operators we can write all the d — 5 contributions for 
each channel. The results are listed in Appendix B. Note that we did 
not assume any specific SUSY model or any texture for fermion masses. 
Our analysis is quite general. 

If the Baryon and Lepton numbers are not conserved, there are many 
channels for the decay of the proton: 



where i = 1,2, 3, and 

n^(K-,n-,K*-,p-)e1r 

where j = 1, 2, while for K* only j = 1. Note that there are also new 
channels for neutron decay. 

These operators have been studied on and off for the last 20 years 
with culminating conclusion that the minimal SUSY SU(5) is ruled 
out [78]. 

In this paper by Murayama and Pierce, the different constraints on the 
Triplet mass are studied. Using the unification of the gauge couplings 
and the proton decay experimental lower bounds they found inconsis- 
tent limits on this mass. 

To give an idea of the procedure used in reference [78], we can use 
the renormalization group equations for the gauge couplings at one- 
loop (neglecting the Yukawa couplings). 
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Now assuming exact unification we get at one loop level: 

3a^ l {jn z ) — 2a^ 1 (m z ) — a{ x {rn z ) = 

1/12 M T m SU SY\ (a 10 v 

— ( — log 2 log J 6.18 

27T v 5 m z m z > 

Therefore it is possible invert the above equation and determine the 
colored Higgs mass. For numerical calculation, they used the two-loop 
RGEs for the gauge and Yukawa couplings between the SUSY and GUT 
scale. 

Knowing the values of the gauge couplings at the scale m z [53], they 
found that the SU (5) prediction of exact unification agrees with data 
only for colored Higgs masses of: 

M T < 3.6 x 10 15 GeV (6.19) 

The second limit on the Triplet mass is computed using the experimen- 
tal lower bound for the channel p — > K + i*i of 6.7 x 10 32 years [75]. 

Computing the proton lifetime due to the d = 5 contributions using 
the methods of reference [79], they found: 

M T > 7.6 x 10 16 GeV (6.20) 

in order to satisfy the experimental bounds. It was assumed nearly 
degenerate scalars at the weak scale, or order 1 TeV in mass. 

Comparing this equation with equation (6.19), they claim that the min- 
imal SUSY SU(5) theory is excluded by a lot. 

Now in order to consider a different scenario and see the possibilities 
to suppress the d — 5 operators, they considered a second case, the so 
called decoupling scenario [80, 81, 82]. 

In this case the first and second generations of superpartners could be 
heavy without severe fine-tuning because they do not affect the Higgs 



76 



CHAPTER 6. PROTON DECAY IN THE SUPERWORLD 



boson self-energy at one-loop level. 

Since the loop factor goes like mg m~ 2 (when mg <C m~), where mg is 
the gaugino or higgsino mass, while m~ is the slepton or squark mass. 
Therefore we can get a large suppression by making the sfermions of the 
first two generations very heavy. For example if we assume m~ 2 ~ 10 
TeV, we get an extra suppression factor of 1CT 2 to the amplitude. 

Now computing the lower bound for the Triplet mass in the decou- 
pling scenario, they found that: 

M T > 5.7 x 10 16 GeV. (6.21) 

Therefore from this analysis they concluded that the Minimal Super- 
symmetry SU (5) is ruled out. 

However in the last analysis, they did not consider the most general 
scenario, as we mentioned in the last section the d = 5 contribution are 
quite model dependent. 

Murayama and Pierce assume the following in their important anal- 
ysis: 

• They found the different limits on the Triplet mass in the minimal 
SUSY SU(5) model without higher dimensional operators in W5, 
which makes wrong predictions for the fermion masses of the first 
two generations. This is not a realistic model, or we could say 
that it is already ruled out. 

• They computed the proton decay amplitudes in a very specific 
SUSY model, where the mixings between fermions and sfermions 
are known, however as we mentioned before in a general SUSY 
model the situation may be quite different. 

• Assuming exact unification of the gauge couplings, they computed 
the limits on the Triplet mass, however if the non-renormalizable 
contributions to Wh are considered, the bounds change. 
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For these reasons we think that the model is not ruled out, in our opin- 
ion before ruling out the minimal realization of the idea of unification, 
it is better to see how constraint the model using the experimental 
bounds. 

In the next two sections we will point out different aspects to be con- 
sidered in order to satisfy the experimental bounds on proton decay. 

6.5 (S)Fermion masses versus Proton de- 
cay 

In this section we will show how to satisfy all the experimental bounds 
on proton decay. As we mentioned before the prediction of the d = 5 
contributions will depend on the explicit form of fermion and sfermions 
mass matrices [83, 84, 85]. 

Now to start our analysis we must write all the d = 5 operators in 
the physical basis, therefore we must go from the flavour to the physi- 
cal basis as: 



F p -> U F F f (6.22) 

F p - U F F f (6.23) 

where F and F represent the fermion and sfermion fields respectively, 
while Up and Up are unitary matrices. 

In general there is not relation between these matrices, therefore we 
expect that in each d — 5 contribution we will have a combination of 
different mixings between fermions and sfermions, which are unknown. 
Note that we only know the mixings between left-handed fermions, 
Vckm = UW and V x = N^E, where N and E are the matrices which 
rotate vp and tp respectively. 

Now the main point in our analysis is, how to suppress the d = 5 
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operators in order to satisfy the experimental bounds. 

Assuming the decoupling scenario, where the most important contri- 
butions come from the third family of superpartners, and studying all 
the contributions shown in Appendix B, we see that the longevity of 
the proton can be achieved by, say, the following conditions at 1 GeV 
[86]: 

D) 3a » (£t D) 3a » (El Echa ~ 
(JV* E) 3a » (Dl Doha ~ {& E) 3a « 
(E/£ J D) 3a ~ (6.24) 

and 

A = e ab (D T C N) a3 (U T A D) 3b » (6.25) 

where a,b — 1,2. Note that [/ and Uc are the matrices which rotate 
u and -u * respectively. We use the same notation for the rest of the 
fermions and sfermions. 

In the above equations we simply mean that all the terms must be 
small. How small? It is hard to quantify this precisely and, honestly 
speaking, it seems to us a premature task. 

Our aim was to demonstrate that the theory is still consistent with 
data and from the above formulae it is obvious. If (when) proton decay 
is discovered and the decay modes measured, it may be sensible to see 
how small should the above terms be. 

Suffice it to say, that a percent suppression of the results in the minimal 
SUGRA (or super KM) case should be enough. This means that on the 
average each vertex should be suppressed by a factor of 1 /3 or so with 
respect to the minimal supergravity predictions. It is very difficult to 
say more: in fact one could be tempted to estimate that for example 
the combinations on the lefthand sides of the above equations need to 
be at least 10~ 2 the same combinations in super KM. However this is 
generally neither necessary nor sufficient. The fact is, that we have to 
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deal with a nonlinear system, since the total decay in a specified mode 
is proportional to the square of a sum of single diagrams, each of them 
is proportional to a combination of unknown mixings. Some of these 
mixings contribute to different diagrams, and some depend on others, 
so the task of constraining them numerically seems exaggerated in view 
of our complete ignorance of all these parameters. What we can say 
for sure is that if each of the diagrams in Appendix B is suppressed by 
a factor of 1/100 with respect to the minimal supergravity predictions, 
proton decay is not too fast and Minimal Supersymmetric SU(5) is not 
ruled out. 

Notice further that the so called super KM basis, in which the mix- 
ing angles of fermions and sfermions are equal, for example Uc = Uc, 
does not work for the proton decay, since eqs. (6.24) and (6.25) are 
not satisfied. If you believe in super KM, you would conclude that the 
theory is ruled out. It is obvious though, from our work, that this is 
not true in general. 

Notice even further, that all the relations (6.24)-(6.25) do not require 
the extreme minimality conditions: A = B_ = Yjj = Yy and C_ = D_ = 
Yd = Y£. More precisely, one can opt for the improvement of the 
fermion mass relations and still save the proton. 

One could worry that the above constraints for the sfermion and fermion 
mixing matrices could be in contradiction with the experimental bounds 
on the flavour violating low energy processes. Fortunately, this is not 
true. Namely, the same conditions (6.24)-(6.25) suffice to render neutral 
current flavour violation in-offensive (of course, the decoupling is nec- 
essary for this to be true). We studied the processes \x — > cy, b — > S7, 
B — B, K — K, etc. It is easy to see, that the combinations which ap- 
pear in (6.24) are exactly the ones that appear in these flavour changing 
processes. So they automatically take care also of these low-energy ex- 
perimental data. The only flavour changing processes that could get 
sizeable contributions are the ones which involve up type sfermions like 
for example D — D or c — > wy. These are not constrained by (6.24), 
but at the same time are not very much constrained by the low-energy 
experiments, so they do not represent a real issue at this stage. 
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Constraints (6.24) are not unique. One can find other relations be- 
tween sfermion and fermion mixing matrices that make the proton de- 
cay amplitude zero or small. However, typically, these solutions can be 
dangerous for FCNC processes, since they do not automatically cancel 
their contributions. So one has to analyze the FCNC processes case 
by case. At the present day status (or ignorance) of proton decay and 
FCNC experiments we believe that this is premature. 

The analysis in Appendix B has been done with the assumption of 
no left-right sfermion mixing, and gaugino higgsino mixing in the neu- 
tralino and chargino sector. This mixing can be included in a perturba- 
tive way, one can show that, up to two mass insertions, the same con- 
straints (6.24)-(6.25) kill all the contributions to nucleon decay. This is 
enough to increase the nucleon lifetime above the experimental limit, 
since each mixing multiplies the diagram by at least 1/10. 
Actually we could ignore the left-right mixing for sfermion proportional 
to the small ratio Mz/m^ in fact, as long as tan /3 > 10, the LR mixing 
can be safely put even to zero without contradicting the experimental 
constraints on the Higgs mass [31] . 

One can also worry about naturalness [87, 88, 89, 90]. Through the 
large top Yukawa couplings, the formula (1.1) becomes here (i = 1, 2, 3) 
(for large tan (3 there are similar contributions of (s)bottom and (s)tau) 



m 2 h « ml + M. [(&U) i3 (tfU) 3i ml + (^ck(4^) 3j m| 

(6.26) 

where and m^c are left-handed and right-handed squark masses. 

Now, for rriq a rhgC ~ 10 TeV (a = 1, 2) in the decoupling regime, 
large (WU)aa or (C/JC/ c ) a3 would imply a small amount of fine-tuning 
(ss 1%) in (6.26). Hereafter, we accept that. No fine-tuning whatso- 
ever, although appealing, to us seems exaggerated; after all it would 
eliminate large extra dimensions as a solution to the hierarchy problem. 
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Strictly speaking, one could then ask why not simply push rriq 3 and 
rrigC all the way up to 10 TeV and be safe? A sensible point, we wish to 
have as many as possible superpartners below TeV and thus hopefully 
detectable at LHC. In other words, all the gauginos and Higgsinos and 
the third generation of sfermions are assumed to have masses lower or 
equal TeV, we only take 2 « 10 TeV. 

In this case, we need to worry only about the third generation of 
sfermions. We also assume light gauginos and Higgsinos, mg ps 100 
GeV. 

The constraints (6.24) can clearly be satisfied exactly by the sfermion 
mixing matrices at 1 GeV. It is reassuring that the sfermionic sector 
does not break strongly SU(2). This is consistent with the SU(2) in- 
variance of the soft masses, which dominate the total sfermion masses. 

The last constraint, eq. (6.25), can even be satisfied in the approx- 
imation C_ — Yd — Ye, which is true in the minimal renormalizable 
model, but at M GUT) not at 1 GeV. 

The relation C_ — Yd — Ye is however not stable under running. To get 
an idea of how big this contribution is at the electroweak scale, one can 
take the approximation that the Yukawas do not run. In the leading 
order in small Yukawas (except for y t ) one gets 



A ~ y c yrV* 3 V 23 V 32 V 21 [1 - {M z /M GUT f* /w *'\ xT 1/M x?xl /a (6.27) 

where V is the CKM matrix and x\ = ai(M z )/au- There is only one 
non-vanishing diagram (the rest vanishes due to (6.24)) and it is pro- 
portional to Vi3A) : fortunately, this seems to be small enough. On top 
of this, in the amplitude the combination (6.25) gets multiplied with a 
combination of neutralino soft masses and m^, which can be fine- 
tuned to an arbitrary small (or even zero) value. And, of course, we 
must keep in mind that Mr can be large and A and/or C completely 
different than Y v (Y E ). 



-1/33 -3 4/3 
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In other words, at this point the proton decay limits provide infor- 
mation on the properties of sfermions and not on the structure of the 
unified theory. 

6.6 The Higgs Triplet Mass 

The determination of the GUT scale and the masses (M T ) of the heavy 
triplets T and T responsible for d — 5 proton decay is one of the most 
important tasks in order to estimate the proton decay amplitude. 

Specifically, if we allow an arbitrary trilinear couplings of the heavy 
fields in E and use higher dimensional terms as a possible source of 
their masses [91, 92], it will turn out that Mt may go up naturally by 
a factor of thirty, which would increase the proton lifetime by a factor 
of 10 3 . 

Now to start our calculation, consider the non-renormalizable terms 
in the superpotential for the heavy sector (up to terms 1/Mpi): 

Of course, if A ~ 0(1), we ignore higher- dimensional terms. However, 
in supersymmetry A is a Yukawa-type coupling, i.e. self-renormalizable. 
For small A (A Mqut/Mpi), the opposite becomes true and a and b 
terms dominate. In this case, it is a simple exercise to show that 

m 3 = 4m 8 (6.29) 

where 7713 and 171$ are the masses of the weak triplet and color octet in 
E. In the renormalizable case m 3 = m 8 . 

At the one loop level, the RGE's for the gauge couplings are 

x _i, 1 / 5 Asusy , 33, Mgut 2 M GUT \ 

a, (Adz) = Oi T , H — In In — In 

1 V ZJ u 2tt V 2 M z 5 M z 5 M T J 
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o^\Mz) = + 



2tt 



25 Asusy , , M G ut , M GUT ^ 
— m — — h In — — him 

6 M z M z m 3 , 



a^(M z ) = ^ + ^(-4 In 



A^y _ 31n r^ + ln M^ ) (g 3Q) 



Here we take for simplicity M GUT = M X y = superheavy gauge bosons 
masses, while at the one-loop level we could as well take A S usy = M z . 
From (6.30) we obtain 



2vr (3a 2 1 - 2e*3 1 - c^ 1 ) 
2vr (5a^ - Sc^ 1 - 2a 3 x ) 



This gives 



-21n^ + ^ln^ + 61n^ 

M z 5 M z 



81n ^w 61n 



m 3 

(^/rnim^M^) 1 / 3 



(6.31) 



Mt 

MquT 



M°(^) 5/2 
\m 8 J 



IV1 GUT 



Mgut\ 



1/2 



(6.32) 
(6.33) 



2m 8 / 

Since, in the case (6.29) is valid, m 8 ~ M GUT /Mpi, we can also write 

Mgut ~ (M^ T ) M pz (6.34) 

In the above equations the superscript denotes the values in the case 
m 3 = m 8 . From (6.29) we get 



M T = 32M£ M GU t ~ 10M^ T 



(6.35) 



Now, ps 10 16 GeV and it was shown last year [78] that M T > 

7 x 10 16 GeV is sufficiently large to be in accord with the newest data 
on proton decay. On the other hand, since we had 



Ml < 3.6 x 10 15 GeV 



(6.36) 
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from (6.35) we see that m 3 = 4m§ is enough to save the theory, even in 
the minimal SUGRA model. Obviously, an improvement of the mea- 
surement of r p is badly needed. It is noteworthy that in this case the 
usual d — 6 proton decay becomes out of reach: r p (d = 6) > 10 38 yrs. 

As we see there are many reasons to believe that the Minimal Su- 
persymmetric SU(5) model is not ruled out. The proton decay lifetime 
depends of many unknown parameters, therefore at the same time us- 
ing the proton decay constraint we could learn about a sector which 
is orthogonal to grand unification. In other words, the improved mea- 
surements of proton decay will provide information about the nature of 
supersymmetry breaking (i.e., the soft masses) and the fermionic mass 
textures. 

After our analysis different groups have been studied the same crucial 
issue in the minimal supersymmetric SU(5) model. Our predictions 
about how suppress the d = 5 contributions has been confirmed in 
reference [93]. In this case they computed the proton lifetime in real- 
istic scenarios, where the needed non-renormalizable operators in the 
Yukawa and Higgs sector of the theory are considered. They showed 
several numerical examples where the experimental bounds are satis- 
fied. 

Therefore from our analysis it is clear that the minimal realization of 
the ideas of supersymmetric unification is not ruled out by the Proton 
decay experiments. 



Chapter 7 
Conclusions 



In this thesis we studied two important phenomenological issues in the 
context of supersymmetric gauge theories. Invisible Higgs boson decays 
into neutralinos and proton decay. 

We started with the study of the invisible Higgs decays into two neu- 
tralinos in the context of the minimal supersymmetric version of the 
Standard Model, considering new one-loop corrections to the neutral 
Higgs boson couplings to neutralinos in the gaugino limit. We com- 
puted these important quantum corrections and included our results in 
the Fortran code HDECAY [39] to compute the branching ratios for the 
neutral Higgs bosons. We focused on the phenomenologically most in- 
teresting case of a bino-like lightest neutralino x? as lightest supersym- 
metric particle, but our analytical results are valid for a more general 
gaugino-like neutralino, irrespective of whether it is the lightest super- 
symmetric particle. We found that these corrections can completely 
dominate the tree-level contribution to the coupling of the lightest 
CP-even Higgs boson. The corrections to the couplings of heavy CP- 
even Higgs boson are somewhat less significant, but can still amount to 
about a factor of 2. Since the CP-odd Higgs boson cannot couple to 
two identical sfermions the corrections are suppressed in this case. In 
all cases the corrections can be significant only if some sfermion masses 
are considerably smaller than the supersymmetric higgsino mass 
The Higgs couplings receive their potentially largest corrections from 
loops involving third generation quarks and their superpartners. In the 
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latter case the corrections are also quite sensitive to the size of the tri- 
linear soft breaking parameter A t (and Ab, if tan (3 ^> 1). 

The possible impact of the corrections on the invisible width of the 
lightest CP-even Higgs boson is dramatic, it could be enhanced to a 
level that should be easily measurable at future high-energy e + e~ col- 
liders, even if the neutralino is an almost perfect bino. This would open 
a new window for testing the Minimal Supersymmetric Standard Model 
at the quantum level 

Turning to applications of these calculations, our results motivated us 
to investigate the effect of these quantum corrections in the elastic 
neutralino-nucleon cross section. In this case we found [28] that these 
corrections might change the predicted detection rate of Dark Matter 
LSPs by up to a factor of two even for scenarios where the rate is close 
to the sensitivity of the next round of direct Dark Matter detection 
experiments. 

A new analysis has been performed after our results were published, 
where also the effect of the subleading loop corrections has been an- 
alyzed. In reference [45] a different group studied the one-loop cor- 
rections to neutral Higgs bosons decays into neutralinos in the general 
case, where we have the decays H® — > XmXn (i=l ,2,3) . They confirmed 
our results computing the branching ratios for similar values of the pa- 
rameters. 

In the second part of our work we focused on proton decay in the Min- 
imal Supersymmetric SU (5) model. We studied the d — 5 operators 
contributing to the decay of the proton, writing all the possible con- 
tributions for each decay channel (see Appendix B) in a general SUSY 
scenario. We pointed out the major sources of uncertainties in estimat- 
ing the proton decay lifetime, as the ignorance of the masses of the color 
octet and weak triplet supermultiplets in the adjoint Higgs, and the un- 
known mixings between fermions and sfermions. Non-renormalizable 
operators are considered in order to correct the relation between the 
fermion masses. 
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We found that the Higgs triplet mass may goes up naturally by a factor 
of thirty, when we allow for arbitrary trilinear coupling in the presence 
of non-renormalizable operators, therefore in this case it is possible to 
satisfy the proton decay bounds, when the mixings between fermions 
and sfermions are known. This is the case of the minimal SUGRA 
model. Alternatively, if the standard value of is adopted, the con- 
ditions for the sfermion and fermion mass matrices in order to achieve 
the proton longevity were found. Knowing all the aspects mentioned 
before we can conclude that the minimal version of the ideas of Super- 
symmetric Grand Unification is not ruled out as claimed before. We 
could say that if proton decay is found, it could provide indirect in- 
formation about the nature of supersymmetry breaking (i.e., the soft 
masses) and the fermionic mass textures. Thus opening a new way to 
test models for fermion and sfermion masses. We could say that our 
analysis is also valid for any Supersymmetric version of Grand Unified 
Theories as minimal SUSY 50(10) [33]. 

Our results motivate the analysis in reference [93], where our predictions 
about how to suppress the d = 5 contributions has been confirmed in 
several numerical examples. 

We hope that these results motivate the realization of new studies in the 
context of supersymmetric gauge theories. The loop corrections com- 
puted by us could be important for future studies in collider physics. 
While our ideas of how to suppress or constrain the d = 5 operators 
contributing for proton decay might be useful for understanding how it 
is possible test the ideas of supersymmetric unification. 
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Appendix A 

Passarino-Veltman Loop 
Formulae 

The explicit form of the Passarino-Veltman loop integrals are: 



C (p 1 ,p 2 ,M 1 ,M 2 ,M 3 ) = (27r/i) 



4— n 



d n l 1 
^D X D 2 D Z 



with 



D 1 =l 2 - Ml + ie 

D 2 = (l+Pif -Mi + ie 

D 3 = (I + Pi + P2) 2 - Ml + ie 

d n l l, 



C lx (p 1 ,p 2 ,M 1 ,M 2 ,M 3 ) = (27T/i) 4 -" / - 

J n 



where: 



m 2 D 4 D 5 D 6 

D 4 = (l- Pl ) 2 - Ml + it 
D 5 = (l- p 2 f - Ml + ie 

D 6 = l 2 - Ml + ie 

where Pi and Mj are the external momenta and the masses respectively. 
For more details see reference [94]. 



91 



Appendix B 

Proton Decay Diagrams 



In this Appendix we present the complete set of diagrams responsible 
for d=5 nucleon decay in the minimal supersymmetric SU(5) theory. In 
our notation T and T stand for heavy Higgs triplets; T and T denote 
their fermionic partners; w^ 1 stands for winos, h +j0 and ft,_ are light 
Higgsinos and V stand for neutral gauginos. 

Decay modes : 

p-+ {K + ,-K + ,p + ,K* + )vi and n — > p°, rj, u, K°, K*°)ui 



OC (D T AU) U , 23 (tfD) 32 , 31 (N T E*) i3 (E T Q T U) 31 



<x (D T AU) n>21 (N T E*) t3 (E T C T U) 33 (tf L>) 32 ,3i 



oc (D T AU) 13 , 23 (U^D) 32i31 (U T D*) 13 (D T CN) 3i 
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where i = 1,2, 3. 
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